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SYNOPSIS 

'Microprocessor Application in a Fire Control System’. 

A thesis submitted in Partial Fulfilment of the requirements 
for the degree of MASTER OF TECHNOLOGY by lieutenant Commander 
Jagdish Anand I.N., to the Department of Electrical Engineering, 
Indian Institute of Technology, EANPTJR, in August, 1977. 

This thesis describes the Fire Control Problem as 
applicable to anti-aircraft firing from a naval vessel. 

Existing systems use analog computation for the solution of 
this problem. A digitalised Error Predictor, using an INTEL 
8080 Microprocessor, to solve the basic fire control problem, 
is discussed. The various assumptions, approximations and 
corrections made are also mentioned. The mathematics of 
computation of time of flight is outlined and its computation 
is presented in details. 

A greater accuracy in the system is introduced by 
incorporating in the predictor an error correction technique, 
wherein, an error function for range, elevation and bearing 
is continuously computed and applied to the computed future 
co-ordinates to give a more accurate prediction. This enhances 
the accuracy of the system and increases the probability of 
hit. The computation of the error function is made possible 
by storing the predicted future co-ordinates at very short 
intervals of time and then comparing these with the co-ordinates 
of the actual aircraft positions after the elapse of time of 
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flight intervals to get the rror values, ihes e error 
values are put on stack and are used to fit a curve from 
which the error functions are computed. 

A new dimension to the problem is the implementation of 
fire control system using Microprocessors. An introduction 
to Microprocessor architecture, constraints due to their use, 
description of INTEL 8080 are included as part of the thesis. 
The Microprocessor system software is dealt with in details. 
Assembly language of INTEL 8080, with special reference to 
instructions extensively used is discussed. Special algori- 
thms due to the use of Microprocessors for computation of 
square roots and curve fitting are also included. Fortran 
Programs and the 8080 assembly language program are included 
in the Append! oes. 

Results show that errors in prediction are considerably 
reduced and that a Microprocessor can easily be adapted to 
the solution of a real time application. A dedicated Micro- 
processor, supported by a small amount of specialised inter- 
face hardware and a well designed software support, is shown 
to offer a viable alternative to special purpose discrete 
logic systems. The low cost of microprocessors provide a 
significant economic advantage to systems of this type. 
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INTRODUCTION 

The main aim of any weapon control system is to fire a 
shell(s) from the gun(s) that the system controls so as to 
destroy the approaching target. Underlying this situation is 
a certain element of probability : the target position at the 
time of hit must be predicted, and there is a certain interval 
between the launching and the hit-called the time of flight _ 
during which the projectile is wholly under the influence 
of natural phenomenon. The computation of the time of flight 
and the accurate prediction of the target future position, 
in essence, constitute the fire control problem. 

Existing systems of the Navy use analog computers for 
the solution of this problem. The aim of this thesis is to 
suggest a digitalised error predictor, using an INTEL 8080 
Microprocessor as its main computing element to solve the 
fire control problem. The predictor described incorporates 
an error correction technique which enhances the accuracy 
of the system and increases the probability of hit. The 
Microprocessor system software is dealt with in full details 
and the system program in the assembly language of 8080 is 
also presented* 

The thesis is divided into six chapters. Chapter 1 
describes the Fire Control Problem. It deals with the 
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assumptions pertaining to the movement of the target, the 
motion of the projectile, choice of the co-ordinate system 
and outlines a method for the computation of time of flight. 
Chapter 2 goes into the details of the time of flight"' com- 
putations. ft discusses the corrections that must he applied 
to- the standard trajectory solution and outlines the mathe- 
matics of time of flight computations. The time of flight 
routine used in the system is discussed. 

Chapter 3 discusses the improved prediction technique. 

It outlines the general representation of the system and 
lists the inputs to and outputs from the predictor. It 
describes the computation of the prediction error and the 
various techniques to correct this error. A detailed descri- 
ption of the error correction technique used is presented. 
Chapter 4 deals with the system design. It presents an 
introduction to general Microprocessor architecture, micro- 
processor imposed constraints, the choice of a microprocessor 
for the system, and a brief description of Intel 8080 micro- 
processor. It discusses the system requirements and outlines 
the proposed system using 8080 and other Hardware modules. 

Chapter 5 is a detailed chapter on System Software. In 
this the 8080 assembly language with particular reference to 
instructions extensively used in programming the system are 
explained. The special algorithms used, the memory 
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requirement and timing constraints are discussed. A system 
flow chart is presented and explained. 

Chapter 6 concludes the thesis and outlines the scope 
for further work. The Fortran Program with sample results 
and the 8080 Assembly language Programme are listed in the 
Appendices. 
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1. FIRE CONTROL PROBLEM 

1.1 General Introduction 

Eire Control Problem involves the launching of a 
projectile from a weapon station (that may be moving) at 
a target (that may also be moving) so as to score a hit on 
the target. To achieve this, the first object in all cases 
of gunnery against a moving target is to fire a projectile 
in a direction and at a time to cause the projectile and the 
target to reach a certain point simultaneously. This point 
is called the future position of the target and will be 
denoted by P. The position of the target at the moment of 
firing is called the present position and will be denoted 
by P, The calculations necessary to achieve this object 
constitute the prediction problem in its widest sense, and 
the computing device that sonves this prediction problem is 

called a PREDICTOR. 

0 

A n Important function for any prediction system is that 
of computing the time of flight t^ , the interval between the 
launching of the projectile and the hit, during which the 
projectile is wholly under the influence of natural phenome- 
non. To compute this accurately it is necessaiy to have 
complete knowledge as to how both the target and the proje- 
ctile will move during this interval. Because complete 
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information can nefer Toe available (and if Jit were, the 
necessary computations would be so complicated as to be 
impracticable), it is necessaiy to make cejftain assumptions 
about the movement of both the target and the projectile. 

1.2 Assumptions Pertaining to Movement of Target 

The predicted motion of target is based on measurement 
of its motion prior to the moment of firing and on assumptions 
as to how it will move during the time of flight. The problem 
in this case is how to find the spherical polar co-ordinates 
(B^, S f ,) of the future position P, allowing for target 

movement during the time of flight t f . The information 
available consists of the spherical polar co-ordinates 
( , S , B ) of the present position P and their corresponding 

• e • 

rates of change (R, S, B) . These are assumed to be available 
continuously . 

Systems that were designed till now have made use of 
one of the following set of assumptions as to target movement 
during the time of flight. 

a) The target travels at a constant speed in a straight line 
at a constant height. 

b) The target travels at constant speed in a straight line, 
but not necessarily at a constant height* 

c) The target travels at constant speed, with constant rate 
of change of height and constant rate of change of course. 
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d) The target has oonstant vector acceleration. 

Of these four sets, though the condition in ( a), is the most 
restrictive it is the one that is most commonly fulfilled in 
normal flight. It can usually be expected to hold except 
when the target is taking avoiding action of some kind or 
mano euvvmng for attack. 

Assumption (b) is the most general for which the target 
vector velocity is constant. It is the most common assumption 
in use as the system is operated with a radar and hence Rp, R, 
Bp, B, Sp, § are continuously available. Assumptions (c) and 
(d) permit a curved target track, and although more accurate, 
involve measurement of second derivatives of suitable co- 
ordinates which are not accurately measurable in the present 
existing systems in the country. 

Bor the design considered, assumption (b), i.e. the 
target is flying in a straight line at constant speed, but not 
necessarily at a constant height, is assumed. 

The success of the prediction will depend very much on 
the fulfilment of this assumption by the target. A pilot 
who goes into or straightens out of a turn or dive during 
the time of flight, whether this is done by chance or is 
deliberate avoiding* action, can always reduce very considera- 
bly the chance of a hit, however elaborate the prediction 


formulae 
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1.3 Motion of the Projectile 

The predicted movement of the projectile is based partly 
on experiment and partly on theory, and the accuracy of this 
prediction is naturally dependent on the reliability of both 
of these. It is even more dependent however on the reprodu- 
ction in practice of the conditions under which the experiments 
are conducted or on which the theory is based. These condi- 
tions are liable to variation in many details, such as the 
state of the atmosphere, the speed and direction of the wind 

•pi 

at various heights, wear of the gun at the instant of firing, 
abnormality of the projectile etc. It is possible to allow for 
some of these, but their more or less infinite variety renders 
complete correction impracticable. 

The most important effect for which allowance must be 
made in prediction is the target movement. Thus, this 
problem has been considered separately and corrections due 
to ballistic abnormalities are added on to it subsequently. 

The computation of the time of flight therefore incorporates 
standard conditions in vacuo with standard muzzle velocity 
of the gun. The corrections applied are discussed under a 
separate heading. 

1.4 Choice of Go-ordinate System 

The three basic sets of co-ordinates are : 

(a) Spherical polar co-ordinates - slant range (R), bearing (B) 
and angle of sight (S) 



8 


(b) Cylindrical polar co-ordinates: plan range (r), bearing 
(B) and height (H) 

(c) Cartesian co-ordinates: eaatings (E), northings (N) and 
height (H). 

Referring to Figure 1.1, it is seen that thqy are conne- 
cted by the formulae, 

r = R Cos S 

H = R Sin S 

E = r Sin B = R Cos S Sin B 

N = r Cos B = R Cos S Cos B 

Eire control systems are distinguished by the particular set 
of co-ordinates that are used for the main prediction compn- 
tations. The system designed uses the spherical polar co- 
ordinates as spherical polars are the quantities that are 
measured by a conventional director system. Further the two 

angles required for aiming a gun correctly (B^ anddb) are in 

1 t not 

effect spherical polar co-ordinates, so if these have/been used 
for prediction it will be necessary to convert back. 

1.5 Computation of Time of Flight 

In order to predict the future position the time of 
flight must be known, but the time of flight cannot be 
calculated till the future position has been found. The 

problem is diagramatically represented in Figure 1.2 and 

' ' . ' # " ' ' • . 

can be represented as a set of two equations. 










( 1 ) 

( 2 ) 
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F = f ± (P, t f ) 
t f = f 2 (F) 


The first of these is a kinematical equation determined by 

the target. It gives the position of the target after a 

time t^ and is of course in very much simplified form here. 

It should in fact be two equations, to give and and 

it involves not only the co-ordinates of P but also their 

derivatives. Equation (2) is based on the motion on the 

motion of the shell and is therefore ballistic. These two 

equations can not be explicitly solved and thus an iterative 

method is used to solve these. A trial value of t„ = t 

f p 

(the time of flight to the present position P, approximated 
to present range/muzzle Velocity), is selected. This value 
of t is used to obtain a future position F by means of 
equation (l). The co-ordinates (R^, S f ) of this position are 
substituted in equation (2) to obtain a new trial value t^. 
This cycle of operation is repeated until no appreciable 
changes are obtained in the values of and t^, which are 
then the required values. It can be proved that this method 
is always convergent, and for the cases studied, the algorithm 
used has been such that the final values are obtained in 10 
iterations at start and ia ^7 Iterations for each subse- 
quent values at intervals of 0.1 or 0.2 second. 
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2. DETAILED TIME OF FLIGHT COMPUTATION 

2.1. M athematics of Computation of. Time of Flight and Tangent 
Elevation 

In the Figure 2.1, GT is the initial direction of motion 
of shell and TF is a vertical from T (drop height or gravity 
d'rop) . 

Thus the future position F can he defined by means of the 
lengths GT and TF. For a general trajectory these are deter- 
mined from two differential equations. Owing to the compli- 
cated form of the air resistance law these equations can not 
be solved explicitly. If however the air resistance is 
neglected or is taken to be proportional to the velocity, 
explicit solutions can be obtained and although these are 
not accurate enough to be used as approximations in the 
general trajectory, they do provide suggestions as to the forms 
that emperical formulae might take [l] . For the system under 
consideration, for simplicity, the trajectory is considered 

in vacuo, with no air resistance. Under these conditions 

2 

GT = Vt, and TF = -J- gt , where V is the initial velocity 
(muzzle velocity). 

From Triangle GFT, we have 

Sin E = H Cos S f = || Cos S f 
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Figure 2.1 

General Trajectory Representation 
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.'vs tiie tangent elevation B is small, it is equal to 
f 't f ) Cos The function f(t~) is of course the tangent 

el ovation B q for zero angle of sight and time of flight t^. 
Thus , 

SIN 1 = E = B Cos 

o i 

In the case where the air resistance is pz’o portions! to 
velocity, say W per unit mass where V is the instantaneous 
velocity and fc is a constant, it can be proved that 

GT = £ (1 - e~ kt ) 

IE = ( kt _ i + e “ kt ) 

It 

and similarly applying the sine formula, for small values of 
B we have 

E = E Oos S„ 
o f 

Although the justification of this approximation is valid 
only in the case of vacuo or when air resistance is pro- 
portional to velocity, it is found that a similar formula 
can be used with a general resistance law [l], A graph 
showing the tangent elevation for different values of 
of an existing gun is shown in figure 2.2. The dotted 

curve shows values of E Oos S_. It can thus be seen that 

o ± 

this approximation is reasonably justified as the variations 
from actual is nominal. 



Figure 2.2 

Graph of Tangent Elevation with Fang 
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Prom Figure 2.1, we have. 


h = ® - SI - s f + 3 > 


Hew Op s gf 
cTos" 3- 


Oos ( 3^ + Si) 


Oos b_ 


Oos E — Tail S^ Sin E 


As S' is small this can be approximated to be 


Oos d 

Oos = 1 - Tan E q Oos S^ 

= 1 - E‘ Sin S- 
o x 

Therefore, R^ = G-T (1 - E q Sin S^) 

= Vt (1 - ^ Sin S f ) 
i.e. 2 R^ - 2 7^ - gt 2 Sin S f 

Therefore, g Sin t 2 - 2 Vt + 2 R^ = 0 (l) 

This gives a quadratic in t and corresponds to the second 
equations in the time of flight computation as depicted in 
the previous chapter. Restricting to first order approxima- 
tions the first equation is 

R f = Rp + *R t f (2) 

Prom equations (l) and (2) we iteratively get the time of 
flight and future co-ordinate;, of the target. The future 
co-ordinates for elevation and beering are obtained by 
computing the vertical and azimuth deflection. 



(D v ) and (D a ) respectively. These are given by the 
relation [l] 
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Sin(D v ) 

Sin(D a ) 


i S t 

Ep t f 

Ts J_ 




!os S 


P 


Oos 



S f “ S p + D v 
' and B f = B p + D & 

These co-ordinates are then corrected for the following 


(a) Corrections for non-standard ballistic conditions 

(b) Drift corrections 

(c) Corrections due to the effect of wind on shell 

(d) lateral and vertical convergence due to the physical 
displacement between the directing sy stars, and the guns 
(Refer Appendix A). - 


These corrected co-ordinates are then fed to the guns 
as 'Gun orders' and the gun is therefore 'trained * and 
'layed' to fire at the future position of tile target. 


2. 2 Correction for non-standard ballistic Conditions 

These can be separated into internal and external 
effects. The former include wear of the bore and variation 
of charge temperature, whose effect is to change the muzzle 
velocity. The latter arise oil account of variation of 
external temperature, humidity and barometric pressure and 
. cause a change in ballistic coefficient. 
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Referring to figure 2,1 and equations in last section 

it is seen that a change in the muzzle velocity Y produces 
a proportional change in G-T , but no change in If. This is 
exactly true in case of motion in vacuo and with resistance 
proportional to the velocity and is approximately true 
for the general trajectory. The change is represented as a 
'percentage = 6V. The resistance of the air to the motion of 
the shell depends on the velocity relative to the air and also 
on the shape and calibre of the shell and the density of the 
air. It is in fact inverse^ proportional to the ballistic 
coefficient, defined as 

Mf /KcT d 2 

where K = coefficient of shape 

o' = coefficient of steadiness 
d = calibre 
M = Mass of shell 

f = a factor which allows for the variation of the 
air density with height (unity for zero height) 

T = coefficient of tenuity, the density of the air at 
sea level. It is a function of height. T depends on the 
barometric pressure, absolute temp., and humidity at sea 
level and can be calculated when these are known. The basic 
trajectory is calculated [l,2] for standard values of 
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pressure, temperature and humidity, and variations from these 
values when set in the predictor are combined to give a change 
60 in 0 (expressed as a percentage). Limiting the corrections 
to first order approximations, the corrections to time of flight 
and elevation are linear functions of 6 V and 60 i.e. 

St = A5V + B6C 

and 6S f = A’SV + B’60 

where A, B, A 1 , B* are functions of the co-ordinates of F. 

Since these corrections are small, relatively crude emperical 
approximations can be made. For the system under design these 
were simulated using random number generation procedures. 

2.3 Drift Correction 

This is due to the spinning action given to a shell due 
to the rifling in the muzzle of the gun. This causes a pre- 
cession in the shell making the axis of the shell point 
away from the direction of motion. Drift is usually repre- 
sented as an angular displacement, and is corrected by applying 
an equal and opposite deflection in azimuth or in lateral plane. 
The azimuth deflection for drift is a function of the time 
of flight t^ but is independant of </. The function f(t^) 
is not in fact linear, but since drift is a small angle, a 
linear approximation is usually adequate. The azimuth deflection 
for drift is approximately proportional to tangent elevation 
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at zero angle of sight [l]. The corresponding deflection 
in the future lateral sight-plane is therefore given by 

Drift = Drift constant x B 0 os S„ = Drift constant x 

o i 

Tangent hlevation 

. 4 Corrections due to Wind 

In considering the wind effects it is assumed that tho 
wind is constant and horizontal. The direction from which 
it is blowing is usually represented by the angle WI (called 
the Wind Inclination) measured clockwise from the vertical 
plane of fire. 


If W = Wind speed 

Wr = W Cos (WI) is the component along the vertical 
plane of fire 

WI = W SIE (WI) is the component across the vertical 
plane of fire. 

It is found that due to the wind the appropriate corrections 
are [lj . 

(i) 6 R = Wr t^ Cos (to range) 

(ii) &<p = - (Wr t_~ Sin<fc»)/R.£ (to elevation) 

The range correction is affected by corrections (-■t/$R) 6 R 
and (^VyCR)^fi. to time of flight and elevation respectively. 
These work out to dt^, = Wr t^ Cos S^/V and 6 dv - 
-(Wr Cos S^/V Sin cj> Coscp ) . 
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Like the other corrections these are also relatively small, 
so that the various functions of P can be represented by 
emperical approximations, such as linear functions of range 
or time of flight, obtained usually by graphical methods. 

In considering wind effects, we have assumed constant 
wind speed and direction at all points of the trajectory. 
Through sufficient metrological information an 'equivalent 
constant wind 1 (defined as the constant wind which will 
produce the same displacement in magnitude and direction as 
actual air motion) is determined and used in the computations. 
In all the cases considered wind speed considered is the 
equivalent constant wind. 

2.5 Tine of Plight Routine 

1 flow chart for the computation of tine of flight is 
shown. in Figure 2 . 3 . This TIMS routine is called by the main 
Prediction Error Program to compute the tine of flight and 
the future point of impact co-ordinates of the shell and the 
aircraft, if the firing was to take place at the present 
instant. The tine of flight computation is done iteratively 
within the allowable tolerance of 0.05 seconds in tine and 1 
foot in range. Por the cases studied this routine is found 
to converge in 10 iterations, u safety margin is provided 
by allowing the routine to run to a maxinun of 15 iterations. 
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Time of Flight Flow-Chart 
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The time of flight is dependent upon the future range 
and angle of elevation. The elevation at which a gun can 
fire is generally limited between 10° and 80° as otherwise 
the ship's super-structures come- in the way. To prevent a 
gun from firing outside the limits, mechanical switches 
which interrupt the firing circuit are incorporated in the 
guns. In the iterative calculations of time of flight an 
attempt is made to limit the gun elevation within 10* to 80° 
when the target is in the firing range. 

In the design of the systan there exists a time delay 
in determining the correction function co efficnets. This 
time delay is equal to the initial time of flight, when 
the target is first detected. To minimise this delay the 
correction function coefficients should be obtained earlier. 
This is achieved by assuming a hypothetical muzzle velocity 
which is much higher than the actual muzzle velocity. A 
muzzle velocity of ten times the muzzle velocity is used for 
this purpose. This is however, done when the target is 
detected at a long range (^>.20,000 yards), as it is then 
that the delay is large. The use of the correction coeffi- 
cients determined with the hypothetical shell velocity, 
marginally effect the system performance [ 3 ]. Thus these 
values are used till the correction coefficients can be 
determined with actual muzzle velocity. When the range is 
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-^15,000 yards the actual nuzzle velocity is used for all 
computations. Thus within the firing range of the gun the 
correction coefficients computed with the actual muzzle 
velocity are always used. 

The corrections discussed in the preceeding sections 
have been included in this routine. In the microprocessor 
implementation of this routine, a special algorithm to 
determine square roots in solving the time of flight qua- 
dratic is incorporated. This is discussed in detail in 
Chapter 5» Section 5. 3.1» 



4 



5. IMPROVED PREDICTION TECHNIQUES 
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A simplified fire control system is shown in Figure 3*1. 

It is assumed that the system is stabilised i.e. the gun and 
the radar antenna are mounted on stabilised platforms and all 
data available is stabilised. This is achieved by applying 
correcting signals from the ship’s gyro stabilisers. It is 
also assumed that the gun and the Radar antenna are theoreti- 
cally at the same point and plane. This is achieved by applying 
Dip, Lateral and Vertical Convergence corrections to the angle 
at which the projectile is to be fired. ( These terms are 
explained in appendix A.) 


3.2 Predictor 

Out of the four sub-systems constituting the fire control 


system, the one which is of interest to us is the Predictor. 
The inputs available to the Predictor are : (A) Through the 
Radar and TRACKER (device that generates the necessary signals 
for automatically tracking the target, once the radar is 
’locked-on', to the target). These are : 

1. Present range of the target (R^) 

2. Present bearing of the target (B^) 

3. Present elevation of the target (Sp) 

v * 

4. Rate of change of range of the target (R) 

5. Rate of change of bearing of the target (B) 

Rate of change of elevation of the target (S) 


6 
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(B) The set of the variable inputs that is manually fed to 

the predictor. These are : 

1. Correction for change in Muzzle velocity (6V) - a 
factor determined from manuals particular to the gun 
and dependent on the number of rounds fired since the 
last calibration to the standard muzzle velocity [1]. 

2. Correction for change in Ballistic conditions (6C) - 
dependent upon atmospheric conditions. 

3. Equivalent wind velocity, W 

4. Wind Indination, WI. 

(C) Set of fixed inputs for gun orders. These are : 

1. Height difference between gun and director (h) - 
to compute DIP. 

2. Horizontal distance between gun and director (d) - 
to compute lateral and vertical convergence. 

The outputs from the predictor consist of : 

1. Gun training (B^) 

2. Gun Elevation (S ) 

O 

3 . Time of Plight of Shell (t^) - this is necessary for 
systems that employ time -set fuzes for shells, to 


set the fuze number 
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These inputs and outputs are diagrama tic ally shown in Fig. 
Figure 3 . 2. 

The function of the predictor is to compute B , S and 

S § 

t f from the available inputs. The equations used in the 
computations make use of some constants, which in real life 
situations are slowly vaiying variables. For simplifying 
the predictor design, these are assumed to be constants. 
Emperical relations derived after extensive trials are also 
used. Such imperfections introduce error in the system and 
reduce the accuracy of the system. The error can be divided 
into two parts ; 

(1) Error in prediction of future position of aircraft. 

(2) Error in prediction of shell trajectory. 

The latter error can be reduced by collecting sufficient 
data regarding shell trajectories through a number of experi- 
ments and trials under vaiying circumstances, and deducing 
emperical relations that can be applied as correcting functions 
at the time of firing. A correction technique to reduce the 
error in prediction of future position of aircraft is explained 
in the next section. 

3*3 Correction Technique and Methods 

The predicted future positions of the aircraft at regular 
time intervals is stored to find the error in prediction at 
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various time intervals. The error in prediction is determined 
by comparing the actual aircraft position with that of the 
predicted value after the elapse of the time of flight inter- 


val. If the error in prediction of the range is DR^ when 
the predicted future range is R^ for various i, then DR^. 
is expressed as a function of R^. From a set of values of 




a curve fit routine is used to compute the coefficients 


of this function. Once the coefficients are known, for any 


predicted future range R^, the error in prediction DR^ is 
computed. This is applied as a correction to R^ to reduce 
the error in prediction of - R fx * Similar logic applies to 
correction in elevation (DS ) and bearing (DB ). 

.X JC 


Depending upon the behaviour pattern of the error values 
with respect to range , any of the following methods can be 
used as the curve fit routine to obtain the correction fun- 
ctions : 


1. No correction ; If the error values are always within 
the required tolerance, no correction i« necessary. This 
condition is well nigh impossible. 

2. Constant correction ; If the error values are more or 
less constant (within tolerable limits) over a length 
of time, the predicted values are corrected by this 
constant to give the corrected predicted co-ordinates. 
This is the simplest method to implement in a micropro- 
cessor. However, for the cases studied this has been 



29 


found unsuitable as although the values are considerably 
constant over adjoining sampling intervals, the values 
drastically differ after a length of time and hence it 
is not suited for continuous tracking. 

(c) Average correction : The average of a group of error 

values is used as the correction function. This is better 
than the constant correction method as the values change 
progressively and this method can be used for continuous 
tracking. This is also simple to implement in a micropro- 
cessor. Although, the errors by this method are brought 
within tolerances in a large number of cases, the method 

fails to reduce the errors in some cases. This method is 

between 

suitable where the variation of error values^successive 
sampling intervals is small, or where the error values 
are large. Its ease of implementation in a micropro- 
crssor as compared to the curve fit method makes it a 
contender for consideration in a microprocessor 
system. 

(d) Polynomial curve fitting : This is used in all cases 

where the variation of error with range is smooth, i.e. 
it can be represented by the form 

2 

DX = a Q + a-^H^ -t* 4* ••••••••«* 

If the variation is linear a first order polynomial 
DX = a Q + a 1 Rp need only be considered), 
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This can however, be generalised to higher order polynomials. 
For the large number of cases studied the variation in error 
has been found to be mostly linear and the error in all 
cases is reduced to acceptable values. This thus, provides a 
more accurate technique . However, its implementation 
in a microprocessor system requires much greater memory and 
this takes more time for execution. This has been dealt with 
in details in Chapter 5 under the section on special algorithms, 

(e) Dynamic prediction (adaptive) : where the error 

DX = f( range) where f is a function that 

minimizes some specified objective function as least-square 

a 

error, f is in general ^difficult function to determine and 
it may not beworth wasting the memory and execution time, 

3.4 Prediction Error Routine 

A flow chart indicating the computation and storing 
of Prediction Error is shown in Figure 3 . 3 . When the predictor 
is switched OH or a fresh target is engaged, the system is 
initialised and the real time *C LOCK is set to zero. It is 
assumed that the manual settings depending on the prevailing 
conditions are made prior to switching. The inputs from the 
radar and tracker are sampled and the predictor computes the 
future position of the aircraft and the time of flight of the 
shell incorporating all corrections as discussed in Chapter 2 
(Time routine). The computed predicted co-ordinates are 


Figure 3 * 3 

Predictor ,3rror Flow Chart 
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stored as PRF(l), PSF(l) , PBP(l). The time of flight tfj. is 
added to the clock and stored as TOF(l). Now after the 
sampling interval dt , the inputs are sampled again and the 
future position and time of flight of shell, t^ computed and 
stored as PBF(2), PBF(2), PSF(2). The CLOCK is added to t f2 
and stored as T0F(2). This is repeated for every sampling 
time of 0,1 sec or 0.2 sec. 

When the real time CLOCK equals T0F(1) the radar and 
tracker inputs are sampled to find the actual position of the 
aircraft as given by R^, and B^. These are now compared 
(in the ERBOR routine) with PRF(l) , PSF(l), PBF(l) respectively 
to get the errors DR(l), DE(l) and DB(l). Again the inputs 
are sampled when the CLOCK equals T0F(2) to get the errors 
DR(2), DE(2) and DB(2) and this process continues. 

A stack is set up for DR, DE, DB and the RANGE (R^) for 
which these values are determined. A curve fit routine is 
entered which fits curves between the following pair of vari- 
ables 

1) DR vs R^ 

2) DE vs Rj 

3) DB vs B f 

The coefficients obtained therefrom are used to obtain 
the correction functions for range, elevation and bearing. 
Henceforth, for each input sampling, the future co-ordinates 



33 

PfiF(N), PSF(h) and PBF(N) as well as their corrected values 
PRFO(N) , PSFO(h) and P BP 0(h) are computed and stored. If the 
aircraft is within the firing range of the gun, the firing can 
commence with the corrected values of the future position i.e. 
PKFC(h), P3FC(N) and PBPO(N). 

Sample results of the cases studied indicating the 
variation of error in range, elevation and Bearing without 
correction and with correction are given in Table 3.1(a) to (e). 

The cases chosen correspond to : 

a) Short range, fast approaching target . 

b) Short range, fast receeding target . 

c) Short range, slow approaching target. 

d) long range, fast approaching target. 

e) long range, fast receeding target. 

Dimension for parameters are : 

Range in feet, range rate in feet per second, elevation in 
degrees, elevation rate in degrees per second, bearing in 
degrees, bearing rate in degrees per second, wind velocity in 
feet per second, wind inclination in degrees, elevation and 
bearing errors in minutes, time in seconds. 



Table 
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Sample Data 1 

Constant Values are ¥ = 74.02063 WI -- 21.29906 
EEATB = -700.00 SHATE = 0.500 BEATS = 0.500 

Initial Values are HP = 43000.00 SP =45.00 BP = 135.00 
Initial time of flight = 16.1822 


TIME 

AFT BE 
DETECTION 

U 

28.3 

0 

HANG-B 

ERROR 

530.32373 

24.36743 

ELEVATION 

BEEOH 

-116.60914 

-4.49531 

BEAEING- 

BEHOfi 

61.82303 

-0.95028 


U 

459.77759 

-113.13882 

69.16962 

31.9 


C 

23.94946 

-2.09795 

0.25070 


U 

393.73047 

-109.51945 

76.86166 

35.6 


C 

21.44385 

-0.94767 

0.93338 


U 

332.40503 

-106.51945 

84.99652 

39.2 

C 

19.12573 

-0.82155 

1.26335 


U 

275.98096 

-103.48161 

93.74483 

42.9 


c 

16.87927 

-0.64222 

1.62036 


Fired at target after 42.9 seconds of detection 
Time of flight = 5.1635 seconds 
U = Uncorrected Error 
0 = Corrected Error 






35 


Table 3.1(b) 

Sample Data 2. 

Constant values are W = 88.70099 ¥1 = 23.89209 

RRATE = 400.00 SRATE = -0.500 BRATE = 0.050 


Initial values are HP = 21000.00 SP = 40.000 BP = 180.00 


Initial 

Time 

of Blight = 

11.8073 


TIME 

AFTER 

DETECTION 

U 

27.0 

G 

RANGE 

ERROR 

747.59446 

5.70386 

ELEVATION* 

ERROR 

-170.73746 

-4.38291 

BEARING 

ERROR 

49.70097 

0.69804 


U 

827.81055 

-185.84440 

46.81563 

33.0 


0 

7.69238 

- 8.20768 

1.01285 


U 

908.39844 

-203.58168 

44.11432 

39.0 


0 

7.66748 

-12.29838 

1.18869 


u 

988.87207 

-224.77701 

41.59079 

44.9 


0 

6.91016 

-18.17242 

1.38430 


u 

1068.72461 

-250.54949 

39.24096 

50.8 


0 

2.54883 

-21.26823 

1.12346 


Fired at target after 50.8 seconds of detection. 
Time of flight = 22.4108 seconds. 







Table 3.1(c) 
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Sample Data 3 

Constant Values are W = 89. 6 3019 tfl = 28.03121 

RRATS = -200.00 SRATE = 1.000 BRATE = 1.000 

Initial values are HP = 18000.000 SP = 45. 000 BP = 260.00 
Initial time of flight = 7.5805 

TIME AFTER RANGE ELEVATION ' BEARING 

DETECTION ERROR ERROR ERROR 


U 292.854 86 -I36. 33244 103.85501 

14.9 

c 2.33337 0.00430 5.03087 


19.4 

U 

c 

240.74109 

5.78442 

-129.66821 

-3.11157 

126.80995 

11.14218 


u 

192.00696 

-124.49663 

158.66980 

23.8 


c 

5.79224 

-2.52230 

14.91912 


u 

146.97839 

-119.95498 

207.26841 

28.3 


c 

5.32153 

-1.94387 

22.34560 


u 

106.10266 

-116 . 29331 

291.49377 

32.8 


c 

4.50317 

-0.98105 

39.64620 



Table 5 . 1 (d) 
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Sample Data 4 

Tonstant values are W = 88.98053 WI = 61.56277 

RRATE = -600.00 SRATE = 2.000 BRATE = 2.000 
Enitial Values are BP = 90000.000 SP = 25.000 BP = 90.000 
Initial Time of Plight = 3.9071 


DIME IP TEH 
DETECTION 

U 

3.2 

0 

RANGE 

EBR0B 

215.00098 

-0.96289 

ELEVATION 

ERROR 

14.59560 

-0.34180 

BEARING 

ERROR 

-11.43042 

0.21168 


U 

192.73145 

16.09946 

-10.33763 

L3.0 


C 

-2.25098 

-1.01016 

0.69916 


u 

169.17676 

17.08687 

-8.13957 

L7.8 


c 

-0.86523 

-0.35363 

0.90122 


u 

145.15625 

17.74512 

- 3.61605 

22.5 


c 

-0.10449 

-0.14215 

1.46479 


u 

121.45605 

18.16956 

7.30400 

27.3 


c 

0.34570 

-0.21747 

6.38219 


Fired at target after 27.3 seconds of detection 

Time of Plight = 3.2233 seconds (with increased muzzle velocity) 
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Table 3.1(e) 

Sample Data 5. 

Constant Values are ¥ = 88.04145 WI = 58.61721 

RRATB = 500.000 SRATB = -1.500 BRATE = -0.050 
Initial Values are HP = 75000.00 SP = 45.000 BP = 1.000 
Initial time of flight = 3*3225 


TIME AFTER 
DETECTION 

U 

7.2 

C 

RANGE 

ERROR 

98.03516 

-0.37012 

ELEVATION 

ERROR 

14.71017 

1.26967 

BEARING 

ERROR 

13.20212 

0.14789 


U 

112.26367 

11.80898 

12.18540 

12.2 


c 

-1.29980 

-1.05926 

0.39178 


u 

124.68555 

6.47944 

11.49490 

17.2 


c 

-1.14355 

-1.76174 

0.18448 


u 

134.93848 

-5.84812 

11. IO 83 8 

22.2 


c 

-2.39355 

-7.52430 

0.33099 


u 

142.69043 

-61.10007 

10.84027 

27.2 


c 

-2.02441 

-40.62914 

-0.10405 

Fired 

. at target after 27.2 seconds, 


Time 

of Flight 

= 3.9124 seconds 

(•with increased 

muzzle velocity) 
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4. SYSTEM DESIGN 

4.1 Introduction to Microprocessor Architecture [6,7] 

A review of the microprocessor is given in this section. 
A digital computer basically involves a control unit, an 
arithmetic and logical unit, and a randomly addressable 
memory. Instructions and data are stored in the memory. A 
microprocessor is an arithmetic and logical unit plus control 
unit realised on a small number of LSI chips. It provides 
the basic arithmetic and central processing unit of the com- 
puter. The instruction decoding, in microprocessors, may be 
done by random logic realized on the chip, or ah alternative' 
is to build a very elementary instruction set, called micro- 
instructions, whose functions are simply related to the hard- 
ware and then to realize each macro -instruction as a ’sub- 
routine’ made up of these microinstructions. 

B.vte Size 

Most microprocessors have 8 bits for a byte. In most 
of the communication systems, display systems, one does not 
generally require more than 8 bit data and further 8 bits of 
data is a respectable level for A/D, D/A conversion. However, 
for complex arithmetic requirements microprocessors with 
16 bits are available. 
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Registers 

Registers correspond to writeable memory and so are 
especially important in microprocessor applications in which, 
most memory is read only. 3 general uses of registers are 

(i) for arithmetic and logical operations 

(ii) for input/output and memory addressing operation 

(iii)for maintaining control over the microprocessor operation 

itself (stacks, program counter an 4, special purpose 
registers) 

One or more accumulators or general purpose registers are 
normally provided and used both for temporary storage of 
data and also to provide source and destinations for opera- 
tions, performed in the ALU. A large number of microproce- 
ssors also have an Index Register to facilitate indexed 
addressing. This can usually be increment ed/decremen ted 
under program control. Some microprocessors are provided with 
only GPRs and these can be used either as . accumulators or as 
index registers. Application programs are generally organi- 
zed in subroutines and interrupt handling routines, the 
former to modularize the software and the latter to synchro- 
nize the microprocessor with the external inputs. Both lead 
to the need to interrupt normal processing sequence, process 
a routine, and then return to the point of interruption. 

Ror this, microprocessors provide some kind of a stack of 
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registers which, facilitate the saving and restoring of the 

program counter. In some, a fixed set of registers in the 

microprocessor itself provide this facility (aided by 

special instructions for adding an address to the stack and 

deleting it later). In others, only pointers to the stack 

are provided in registers in the microprocessors, with the 

in 

stack itself maintained^writeable memory. 

Microprocessors also incorporate a status register 
which indicate the internal states of machine. 

Most microprocessors use a microprogrammable controller. 
This is typically organised as shown diagrammatically in 
Figure 4.1 below. 



4 different types of data come into it : 

(i) information from instruction register (IR) 

(ii) timing information (</) (a) square wave output of clock 

(b) external timing signals 

(iii) internal states of machine, i.e. information from 
status register (SR) 
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(iv) external Input /output 

The outputs from the controller consist of 

(1) Enable multiplexers - this sets up data paths 

(2) load signals - this along with (l) above is required 
to enable registers to communicate amongst themselves. 

(3) Signal to arithmetic and logic unit (ALU) to set up 
states 

( 4 ) External Input /output signals. 

Bus Organisation 

In a microprocessor system there is one bus in the system 
with n lines, and devices are 'hung on the bus ' . A typical 
organisation is as shown in Figure 4.2. 

An instruction fetch cycle would thus have the following 
sequence : 

PO contains address of instruction 


(i) BUS < PO 

(ii) MA < BUS 

(waiting period - memory fetch) 

(iii) BUS MD 

(iv) IP « BUS 


Assuming two operand instructions, in dealing with any 
instruction the following 5 pieces of information is required. 
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N -lines 



Address 


DATA 


Figure 4» 2 


Microprocessor Bus 
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(i) Instruction address (normally contents of PC) 

(ii) Operand A 

(iii) Operand B 

(iv) Result 

(v) Next instruction address 

In most microprocessor systems operand A is in memory and 

operand B in some GPR (may Be accumulator). Result is also 

put into some GPR and in simple microprocessor system this is 

is 

also the accumulator. The next instruction address Jin the 
Program Counter (the controller using information from 
Register/ALU loads the next address into the Program Counter) . 

Addressing Modes 

The different addressing modes used in a microprocessor 

are 

(a) Direct 
(h) Extended 

(c) Immediate 

(d) Indexed 

(e) Relative 

One major addressing mode that is absent in microprocessors is 
the indirect addressing ? that is the full address specified is 
the address of the address. However, microprocessors do have 
’Register Indirect’ addressing ,where the register or register 
pair specified contains the address. 



Interrupt System 
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The external interrupt is the main scheme by which 
external events and devices are synchronised with - 

internal programs in a microprocessors. External interrupts 
signal the microprocessor that some event has occurred and 
the microprocessor then executes the application program 
associated with that event. The sequence of events in case 
of an interrupt is 

(i) Store all information, i.e. save all the processor 
registers 

(ii) Service the interrupt 

(iii) Clear DONE 

(iv) Restore all processor registers 

(v) Reenable the interrupt, i.e. execute El (enable 

interrupt) 

(vi) Return to interrupted program 
Memories t 

Memories used with microprocessors could be of any 
v : technology depending on the speed and space requirement. 
The various processes are : 

(i) E-MOS memories 

(ii) P-MOS memories (occassionally - advantage is that 
voltage range can be greater) 

(iii) CMOS - slowest but uses least amount of power 

(iv) Bipolar - faster than EMOS by a factor of 10 
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- consumes much more power 

- much larger in area and volume and therefore 
densities are not as great 

(v) EGL - fastest (memories of order of nanoseconds) 

The different types of memories that would be required 
depending on the applications and cost in mind are : 

(f) BGM - read only memories - usually used for storing 

fixed programs. These are masked at the factory 

(ii) FPROM-Field programmable read only memory 

-cost per bit is 2 to 5 tines that of BOM 

- in this diodes are linked 

(iii)EPBDM- Electronically programmable BDM 

- these can be erased by ultra violet radiation 

(iv)EBPBDM-Electronically erasable programmable BOM 

- these can be erased electrically instead of 
ultra violet light 

(V)EAPBOM- Electrically alterable PBOM 

- these can selectively erase 
All the above are non-volatile memories 

(vi) BAM Random access memory 

- memory which can be written or read 


(vii)SSEAM Solid State BAM 
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4.2 Microprocessor Imposed C o nstraints 

There are several important limitations of microproce- 
ssors which must be considered in the design of a micropro- 
crssor based system. Any system that is designed without 
consideration of these limitations will result in a system 
which is, at best marginal in performance. 

(a) Speed : Compared to either analog or discrete logic 

systems, microprocessors are slow. Gate time for standard 
TT1 discrete circuits are measured in nanoseconds. Even the 
faster E-channel metal oxide semiconductor (NMOS) micropro- 
cessors are much slower. A typical memory to register add 
instruction will require 2 to 10 ^seconds in most micro pro- 
cessors . 

(b) Short Word Length ; Microprocessors typically operate on 
short word lengths. Most microprocessors utilise words 

of from 4 to 16 bits with 8 bits being most common. A word 
length of 8 bits is only sufficient for values upto 256. 

Thus in application, like ours, where angles upto 359.9 
degrees, with resolution to the nearest tenth of a degree, 
and values of range upto 30,000 yards with accuracy of 1 yard, 
multiple word storage of data will be necessary. This will 
require additional time for computation. 
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(c) Limited Arithmetic Capability : Microprocessors have very 

United arithmetic capability. There are currently no 
microprocessors with multiply or divide instructions. Jbr 
applications that are not tine bound one overcomes this 
handicap through software routines to multiply and divide. 
However, in real tine applications, like ours, where the 
execution tine is critical, hardware support in terns of 
multipliers and dividers is essential and is therefore 
suggested in the design. 

4.3 Choice of Microprocessor 

A very debatable point is the preference of a micropro- 

erssor to a discrete logic system. In the systems under design 

for 

a nicroprocessor was found to be essential ^he following 
reasons « 

(a) With a nicroprocessor there is a considerable flexibility 
in the choice of different parameters (like the setting of 
different coefficients). Also as there is flexibility in the 
choice of co-ordinate systems, use of nicroprocessor makes 
things simpler in transforming from any co-ordinate system 

to the desired co-ordinate system. 

(b) The system requires data buffering in that the predi- 
cted co-ordinates are continuously stored and along with 
other changing, parameters are required for continuous 
computation. 
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(c) The microprocessor response is made fast enough with 
the inclusion of hardware modules for time consuming 
floating point arithmetic and trignomotric functions. 

The reasons for favouring Intel 8080 were the availabi- 
lity of product literature, software support, its ease of 
availability in the country, a rich instruction set and the 
fact that it is recognised as an industrial standard. A 
16 bit microprocessors may perhaps have been better suited for 
the arithmetic but working with 4 bytes floating point arithmetic 
on 8080,- with added hardware support is considered adequate to 
meet the timing constraints. An exhaustive evaluation pro- 
cedure was not done because sufficient amount of literature 
about competitive processors was not available. 

4.4 Description of Intel 8080 

The Intel 8080 microprocessor is a 40 pin singLe-chip 
NMOS GPU which contains instruction decoding hardware on 
the same chip as the arithmetic and logic unit and is upward 
compatible with the Intel 8008 microprocessor. The 8080 has 
a 16-bit address bus, a 8-bit bidirectional data bus and 
fully decoded, TTL- compatible control outputs. It can 
support upto 64 h bytes of mixed MM and BOM memory, and 
can address upto 256 INPUT /OUTPUT ports. It has a basic 
cycle time of 2 (isec [5]. 
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Architect -ire 

The 8080 CPU consists of the following functional 
units [ 8] 

(i) Register array and address logic 

(ii) Arithmetic and logic unit (ALU) 

(iii) Instruction Register and Control Section 

(iv) Bi-directional, 3 state data bus buffer 

A block diagram of the 8080 CPU is shown in Figure 4.3. 

1. Register Array and Address logic 

The register section consists of a static RAM array 
organised into six 16-bit registers. These are (i) Program 
Counter (PC) which maintains the memory address of the 
current program instruction and increments upon every instru- 
ction fetch, (ii) Stack Pointer (SP) - maintains the address 
of the next stack location in memory, (iii) - (vi) Eight 
8-bit registers arranged in pairs, referred to as B,Cj D,Ej 
H,L| and WZ. The temporary register pair W,Z is not program 
addressable and is only used for the internal execution of 
instructions. The other six general purpose registers can be 
used either as single registers (8-bit) or as register pairs 
(16 bit.) . ■ ' 
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2. Arithmetic and logic Unit 

Arithmetic, logical and. rotate operations are performed 
in the ALU. Registers which. function with the ALU induce 
an 8-bit accumulator, an 8— "bit temporary accumulator (ACT), 
a 5-bit flag register (zero , o ariy , sign, parity and auxi- 
lliary carry) and an 8-bit temporary register (TMP). The ALL 
is fed by the temporary register ( TMP) and the temporary 
accumulator (ACT) and carry flip flop. The result of the 
operation can be transferred, to the internal bus or to 
the accumulator. 

3. Instruction Register ( XlO aja< 3- Control Section 

The IR is a 8— bit regi® ter which transfers information 

from the internal bus to the instruction decode and control 
section. The output of the decoder, conbinea with timing 
signals from the timing and state control section, provides 
the control signals for the memory, ALU, and data bus buffer 
Outputs from the instruction decoder and external control 
signals feed the timing and state control section which gene 
rates the state and cycle timing signals. Inputs into the 
timing section include cloolc inputs an -d <£>2 (^Ue inputs 
of a two-phase external cloolc ) an d four internal control 
inputs READY, INT (Interrupt) 9 RESET, and- HOLD. Six dire ct 
control outputs provided for* external control of memory and ; 
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I/O devices are INTB (Interrupt Enable), HLDA (Hold Acknow- 
ledge), D3IN (Data Bus IN), SYNO, ¥R (Write), and WAIT. 

4. Data Bus Buffer 

The data bus buffer is an 8-bit, bi-directional, tri- 
state buffer used to isolate the internal bus of the processor 
from the external data bus ( D Q - D^).. In the output mode, 
the content of the internal bus is transferred into an 
8-bit latch. This latch drives the data bus output buffers. 
These buffers are switched off during input or non tran e 
operations. In the input mode, data from the external data 
Tbnii® is transferred to the internal bus. 



INSTRUCTION MITCH, MEMOHT RBAD, HWOHr WRITE, STaCK 
HEAD , blAOK WRITE, INTERRUPT, HALT, INPUT and OUTPUT are the 
possible machine cycles within the processor. Ho instruction 
cycle consists of more than five of these cycles. A machine 
cycle in the 8080 consists of three, four or five states, 
where each state is equal to the tine interval of one period 
of the two-phase clock. Typically the clock frequency is 
2 MHz, which results in a state of 0.5 |is duration. The 
five states are referred to as Tl, T2, T 3, 14 and T5. 
machine cycles do not require the full five states, although 
none requires less than three. Hence a machine cycle varies 
from 1.5 to 2.5 micro seconds. The activities of the 
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different states are summarised in Table 4.1. States Tl, T2, T 3 
are common to all machine cycles. If the processor has to 
wait for a response from an I/O device, then the cycle nay 
contain one or more wait states TW. During these first 
three states, data is transferred to or from the processor. 

Table 4.1 8080 STATS DEFINITIONS 


State 


Associated Activities 


Tl A memory address or I/O device number is placed on 

the address bus, status information is placed on 
the data bus . 

C2 The GPU samples the READY and HOLD inputs and checks 

for Halt instruction 


Tv/ Processor enters WAIT state if READY is low or if 

(optional) Halt instruction has been executed. 


Ts An instruction byte (FETCH cycle), data byte (MEMO BY 

BEAD, STACK READ or INPUT CYCLE ) or interrupt 
instruction ( INTERRUPT cycle) is input to the GPU 
datibiffor a data byte (MMOHY WEEDS 
STACK WHITE, or OUTPUT cycle) is output onto the 
data bus . 


T4,T5 

(opti- 

onal) 


bates T4 and T5 are available if the execution of a 
articular instruction requires them, if not, the 
PU may skip one or both of them. T4 and T5 are 
nly used for internal processor operations. 
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IESTRUGIIOE SET 

The 8080 Instruction Set consists of 78 instructions. 

The instructions are either one, two or three bytes in 

length. Multiple byte instructions are stored in successive 

memory locations. The first byte gives the OP code, whereas 

the second and third bytes, if required, give the data or 

address. One byte instructions include register-to-register, 

memory reference, arithmetic, logical, Rotate, Return, Push 

Pop and Enable or Disable Interrupt Instructions. Two-byte 

instructions include immediate node and I/O instructions. 

Three byte instructions include Jump, Gall, Direct, load and 
ins t ru c t io ns . 

Store/ The instruction set can be broadly divided into four 
groups : 

(a) Data Movement Group - This includes movement of data 
between registers, register to memory and memory to 
register. Eleven such instructions are available. 

(b) Data Manipulation Group - This includes the instructions 
for performing arithmetic and logical operations on data. 
Twenty nine such instructions are there. 

(c) Decision and control group - This group has twenty nine 
instructions which include branch instructions, uncon- 
ditional and conditional (permitted for eight conditions) 
and for subroutines. 
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(d) Stack, I/O and Machine Control Group - Nine instructions 
for performing stack and I/O operations. Instructions 
for pushing and popping register pairs and status words 
to and from the stack are available. 

Addressing Modes ; 

8080 has the following addressing modes ; 

(i) direct addressing 

(ii) register addressing 

(iii) register indirect 

(iv) immediate addressing 

It has no relative addressing. 

4 • 5 System Requirements 

A practical system should have a detection range 

from zero degrees to 559.999... degrees (that is 0 to 

6.2854 radians) in bearing and distance to atleast 

50,000 yards. Assuming an accuracy of upto one minute in 

angular values and of one foot in distances, the system 

-4 5 

should be capable of representing values from 10 to 10 . 
Intermediate calculations could approximately range from 
10~^ to 10 + ^. The ratio of the largest to the smallest 

Ip 

number is about 10 . If fixed point arithmetic is used 

approximately 4l bits words would be required. Since the 
range of numbers is large, floating point representation of 
numbers is indicated. A 6 digit decimal number would require 



57 


approximately 20 bits for a binary representation. Thus two 
bytes of 8080 would be insufficient to contain the number and 
hence 3 bytes (i.e. 24 bits) are used as the mantissa (nor- 
malised floating point). A 6 bit exponent is capable of 

+63 

handling numbers in the range of 2T- ' which is more than 
adequate to handle the range of numbers encountered during all 
computations inthe system. The arrangement adopted for 
representation of floating point numbers in the system there- 
fore is to use 4 contigous bytes. The first byte is the 
exponent and the succeeding three bytes are the mantissa. 

The first bit of the first byte contains the sign of the number 
and the second bit the sign of the exponent. The mantissa is 
normalised and arranged that the most significant bits occupy 
the second byte with the most signif icant bit as the first 
bit of the byte and the least significant bit as the last bit 
of the fourth byte. Examples of 3 numbers in this represen- 
tation are shown in Figure 4.4. 

As no multiply, divide instructions are provided in 
8080, software routines to do this are necessary. Software 
routines to carry out four byte floating point multiplication 
division, addition and subtraction take too much time in view 
of the large number of times these are required in the system. 
The system being time bound to the execution of the entire 
program in less than 100 milliseconds, hardware imple- 

mentation of the above arithmetic functions is recommended. 



(c) 

Figure 4,4 


Floating Point Representation of 3 Sample 
Fumbers 

(a) + 10000 

(b) - 32.125 

(c) - 0.20 
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This should have an execution tine of less than 50 micro- 
second. i ? or similar reasons hardware realisation of the tri- 
gonometric functions (sine, cosine and arc sine) is also 
r o co nn ended. 

6 Proposed System Using 8080 and Other Hardware Module 

block schematic diagram of the proposed fire-control 
system using 8080, is presented in this section. 

In our case, the. three important considerations are as 
follows : 

1. Floating point hardware nodule ; 

Is the fire control problems involves a lot of floating 
point calculations, it is recommended that a floating point 
hardware is included. It interacts with GPU as follows. Since 
a floating points instruction is not a regular 8080 instruction, 
an illegal instruction trap will occur. In the trap handling 
routine, we will check whether it is a floating point- 
instruction or not. If so control is transferred to the floating; 
point hardware. The floating point hardware, is capable : 

of accessing the registers and 'the memory. (This is so 
because at the tine the trap occurs, the execution of GPU ; 

is abandoned) . Once the floating point hardware has comple- ■ j 
ted. the execution, it acknowledges to the control unit of 
GPU of its completion. 



OLOOk 



Figure 4.5 


Block Schematic of the System 
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2. Sin/Oca /ARSIi'T Hardware ; 

We would also require the hardware for finding the SIR/ 

CO S /aRS IH values. Here there is a trade off between the 
execution speed (i.e. series expansion) and the memory space 
available (i.e. table loop-up). In either of the cases, the 
procedure sinilar to the one discussed above could be adopted. 

3. Synchronization 

We have assumed that the range ( R ) , the elevation (S^), 
and the bearing (B ) are available continuously through an 
u/D converter as input in fixed memory locations. However, we 
do not want the values to change during our computational 
cycle. That is we want these memory locations to be locked to 
the input unit during the computation cycle and to be locked 
to the arithmetic unit during the input cycle. If we assume 
the input cycle to be one unit of indivisible instructions, 
then the case is considerably simplified. 

I 11 Figure 4„6 , whenever the execution cycle is 

progressing, the lock is set, so that the manory locations 
cannot be altered. But the A/D convert ion process may go 
simultaneously during the execution cycle. Once the execution 
cycle is completed., the lock is reset and. memory locations are 
altered by the input. The next execution cycle starts, only 
after 100 msec, from the start of the previous execution cycle. 
Once the execution cycle commences, the lock is set. 
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The output of the systen (i.e. tiie gun orders comprising of 
Gun Training, Gun Elevation, and tine of flight for fuze 
setting) goes directly to the gun setting mechanisms. 
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Input - Memory locking arrangement 
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5. SYSTEM SOFTWARE 

5 . 1 Discussion of 8080 Assembly Language with Particular 

R eference to Instruc t ions Extensively Used 

Assembly language programs are written as a sequence 
of instructions which, are converted to executable hexadecimal 
code by a special program called an ASSEMBLER. The statement 
syntax of the 8080 assembly language is as follows : 

Field 1 is, the LABEL field. This is an optional field, 
if 

which present , may be from 1 to 5 Characters long. The first 
character must be a letter of the alphabet or (Wr or? * A 
colon ( :) must follow the last character. 

Field 2 is the 00 DE field. This field contains a code 
which identifies the machine operation (add, jumps, subtract 
etc) to be performed. It contains one of the mnemonics which 
identify each instruction followed by atleast one space. 

Field 3 is the OPERAND field. This provides any address 
or data, information needed by the CODE field. Depending upon 
the code f.ield, the operand field may be absent or may consist 
of one item or two items separated by a comma. The types 
of information that may be requested as items of an operand 
field and the ways in whichit may be specified is listed 
in Table 5.1. 



Table 5-1 

Operand Field Information 
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OPERAND FIELD BEFORMATIOfl - 


Information required 


Ways of specifying 


(a) Register 

(b) Register pair 

(c) Immediate data 

(d) 16-bit menory address 


(1) 

Hexadecimal data 


(2) 

Decimal 

data 


(3) 

Octal 

data 


(4) 

Binary 

data 


(5) 

Location 

counter 

(fO 

(6) 

ASCII constant 


(7) 

Labels as 

signed values 

(8) 

Labels of 

• instructions 

(9) 

Expressions 



Field 4 is the COMMENT field. This is present for the 
programmer's convenience and is ignored by the assembler. 

This must begin with a semicolon ($). 

Instructions Extensively Used 

A large number of the 78 instructions in the Instruct^©, 
Set have not been used. Only 20 of these instructions have 
been used and these are listed and explained below. Many 
instructions from the Data Manipulation Group could not be 
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used as the system uses 4 byte floating point arithmetic. 
Extensive use has been made of the LHLD (load Direct 
Register Pair H,l) and SH1D (Store Direct -Register Pair H,l) 
instructions to manipulate 4 byte numbers. Instructions that 
have been used are 

(A) FOR PRIMARY MEMORY REFERENCE 

1. LHLD Load H and L Direct 

(l) * (address) 

(H) 4 (address + 1) 

2. SHLD store H and 1 Direct 

(address) < (l) 

(address + 1 ) <: (H) 

5 . IDA load Accumulator Direct 
(A) < 4 - (address) 

4. STA Store Accumulator Direct 

(address) (A) 

5. MOV Move Instruction 

The format is MOV dst, Src 

1- - A,B,C,D»E,H,1 or M 

(dst or src not both = M) 



B. FOR MEMORY REGISTER OPERATE 


6. m. Logical Exclusive - or Register or Memory 
with. Accumulator (zero Accumulator) 

This has been used to zero the accumulator since any bit 
exclusive-or ed with itself produces zero 
U) <— (A) ( + ) (A) 

7. GMP Oonpare Register or Memory with Accumulator 
The format of this instruction is 

GMP reg 

The specified byte is compared to the contents of the 
accumulator. The comparison is formed by internally subtractin 
the contents of REG- from the accumulator (leaving both unchang- 
ed) and setting the condition bits according to the result. 

The zero bit is set to 1 if (A) = (r). The carry flag is 
set to 1 if (A) <- (r) 

8. I HR Increment Register or Memory 

Format IITR reg (B,G,D,E,H,L y M or A) 

INR A 

(A) < (A) + 1 

9. DCX Decrement Register Pair 
Format . DGX rp (B,D,H, or SP) 

The 16-bit number held in the specified register pair 


is decremented by one. 
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10. X0HG- Exchange Registers 

(D) v (H) 

(a) < — (i) 

The 16-hits of data held in the H and I registers are 
exchanged with the 16-hits of data held in the D and E 
registers . 

11. SEEL load SP from H and .1 

(Sp) ^ (H,i) 

The 16-hits of data-: held in the H and 1 registers 
replace the contents of the stack pointer SP. The 
contents of H and 1 registers are unchanged. 

12. MV I Move Innediate Data 

Format MVI R, DATA 

(R) -4 DATA (1 byte data) 

13. LXI Load Register Pair Immediate 

Format 1X1 EP, D.iTA (2 hytes data) 

The most significant 8 bits of the 16-bit immediate 
data is loaded into the first register of the register 
pair, while the least significant 8 bits are loaded into 
the next register of the register pair. 

0. FOR BRANCHING- 

14. JMP Jump 

Format JMP ADDR 

(PG) < ADDR 

Jump to instruction with label ADDR 
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. JM Junp if Minus 

Jump if sign lit = 1 

16. JP Junp if positive 
Jump if sign bit = 0 

17 . JNZ Junp if ^ zero 
Jump if zero bit = 0 

18. JZ Jump if zero 
Jump if zero bit * l 

D. FOR stage: OPERATION 

19. PUSH Push. Data onto Stack 
Format PUSH rp ( B, D,H, or PSW) 

The contents of the register pair are saved in two bytes 
of memory indicated by the stack pointer SP. 

(SP-l) < (First register) 

(SP-2) — **— (Second register) 

SP < SP-2 

20. POP Pop Data Off Stack 

Format POP rp (B,D,H» or PSW) 

The contents of the register pair are restored from two 
bytes of memory indicated by the stack pointer SP. 

(Second register) ^ (SP) 

(First register) (SP + l) 


SP 


SP + 2 



INS 'MOTIONS ADDED AS A RESULT OP ADDITIONAL HARDWARE 


PADD Floating Point 4 byte addition 

Format PADD D, ADDR or 
F^DD ADDR, D 

The 32 bit data held in register pairs D and H arranged 
as sign bit plus signed exponent in register D, most 
significant 8 bits of mantissa in register B, next 8 
bits in H and least significant 8 bits of mantissa in 
register L, are added to the address, where ADDR contains 
the exponent and ADDR + 3 the least significant 8 bits 
of the mantissa. This floating point addition is done 
by the additional hardware. The result of the addition 
is in registers D,E,H,L in the form as explained above. 

FSU3 Floating Point 4 byte subtraction 
Format PSUB D, ADDR or 
PSUB ' ADDR, D 

In this case the first operand is treated as the 
minuend (i.e. the number to be subtracted) 

FMUL Floating Point 4 byte multiplication 
Format FMUL D, ADDR or 

FMUL ADDR, D 
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4. S'DIV Floating Point 4 tty tedivision 
Format FDIV D, ADDS or 

FPI7 ADDR, D 

In this case the first operand is treated as the divisor. 

5. 003 Cosine 

Format 003 — 

The Gosine of the floating point number in the 
registers D,B,H,L is returned in the registers D,B,H»L 
in the form as described in 1. 

6. SIR Sine 

Format SIN 

The sine of the floating point number in the registers 
D,E,H,L is returned 'in the registers 

Dj E yH j I 

7. ASS IN Arc Sine 

FORMAT ARSIN 

The Arc Sine of the floating point number in the 
registers D,E,H,L is returned in the registers 

5.2 MACROS 

A macro is a means of specifying to the assembler that 
a symbol (the macro name) appearing in the code field of a 
statement actually stands for a group of instructions. The 
use of macros increases the efficiency of programming and 
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tiie readability of programs. Macros have been used while 
programming the system. Instead of defining these within 
the main program, all the macros used have been defined 
before the start of the main program and these have only been 
referenced or called in the text of the program. 

The format for the Macro Definition is 

Label Code Operand 

name MACRO Plist 

macro body 

SNDM 

’name' is the name of the macro 

’plist' is a list of unquoted character strings identifying 
the dummy parameters appearing in the body of the macro- 
definition. Subsequent macro references or call specifying 
the actual parameters to be substituted for the dummy para- 
meters must adhere to the positionality of the parameters as 
indicated in ’plist 1 . 

5.3 Special Algorithms 

Algorithms used for working with microprocessors should 
be simple, so as to occupy a small portion of memory and 
should take very little time so that the main program is not 
time-bound to the execution of a particular algorithm. In 
the design of this system we have made use of two special 
algorithms that are simple and quick to execute. 



72 


5.3*1 Finding Square Root of a Positive Number 

Finding the square- root of a positive number ii is 

equivalent to finding the positive solution of the equation 
2 

f(x) = X* “4=0. Probably the best-known of all iterative 

methods to solve this is that of Newton Raphson, where you 

generate successive' approximations from the iteration. 

l(x ± ) 

x i+l = x i “ f »(x i ) i#e * for 

the case under consideration 

x. 2 -A 

x i+l = x i " iTi = < x i + £T) /2 

1 1 

The success and rapid convergence of this method is dependent 

upon the choice of the initial guess. If the initial choice 
happens to be near the root the technique converges very 
fast. 

For the system under consideration, this routine is 
required only for solving the quadratic that determines the 
time of flight. 

The quadratic to determine time of flight is 
R f = Vt (1 -§f Sin S f ) 

i.e. g Sin S f 
Therefore, t = 

2 g Sin 


t 2 - 2Vt + 2R f = 0 
2V + 2 - 4 • ( g Sin S f ) . 2E^ 
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The two roots, XL, &2 can then he written as 


XI = TIM + Si4RT/2g Sin S f 
X2 = TEM - 8QRT/2g Sin S f 

where TEM = = — - — ~ 

<i g bm g Sin 

and SQRT = ^4 Y 2 - 8 R f (g Sin S f ) 


When considering the tine of flight we consider the minimum 
positive tine, therefore, if SQRT is less than TERM we consider 
the root X2 if SQRT is greater than or equal to TEM we consider 
the root XI 

if SQRT TEM 

SQfiT = - X2) * 2 S Sin S f 

if SQRT ^ TIM 

Y 

SQRT = (XL - n-oTV-cH* 2 G- Sin S- 
g o.tn £>£•'• x 


X near approximation for the time of flight is simply the 

ratio of the present range to the muzzle velocity, i.e., XI or 

to Rp/Y. 

X2 can he approximate^ Therefore, if M/Y is greater than or 
equal to g — the root corresponds to XL (TERM + SQRT/2g Sin 

S f ) and if Rp/Y is less than the root corresponds to 

X2 (TERM - SQRT/2g Sin S f ). Accordingly the initial guess for 
determining' the root is 



74 


Y ~” w 

- (g"s in -§- “ y ) 2 g Sin , and if 




D. 


P 


g Sin S f * 


X - -X 


Implementing this, ic has been found that upto an 
~6 

accuracy of 10 the square root can be determined in-in 5 
iterations . 


5.3*2 Curve Pit Routine 

Different curve-fit routines for the general purpose 
computers have been available [10 ] . Ihe theory of curve-fit 
routines has also been widely studied in the past. In the 
present context we will develop curve fit algorithm that 
is suitable for implementing in our microprocessor fire 
control system. Instead of going for a more general algorithm 
we attempt to develop a specific algorithm for the reasons 
as stated in the beginning of this chapter. So, instead of 
an elaborate, general curve-fit algorithm, we write a 
specific curve-fit algorithm to suit our specific goals. 

Input ; Input consists of a set of observations namely 
the aircraft range error, or bearing error, or 
elevation error for a set of future range. 

Output t Coefficient of the polynomial, which fits closely 


the calculated errors. 
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From the results of a large number of cases of aircraft 
movements studied, a first order polynomial is sufficient 
to describe the different errors in terms of the future 
range. This can however, be generalised to a polynomial of 
atmost 3rd order if practical observations so warrant it. 


Algorithm : Under least square estimation, we have to 
n . _ 

minimise ^ . (y. - <a. x. ) , where y . r s are errors, and the 
j 1 1 J D 

x. 's are parameters (predicted range at which the error occurs.) 


The minimum of the above function is 0, and this occurs 


when, 


7-5 


£a ± *5 = 0 


Considering a second order polynomial 


y -} = <r j = 1 ’ 2 ’ 3 

J i=0 J 

2 

i . e . 7p — 9 -q "f* a-j x-j -f- a 2 x-^ 

2 

y 2 = a Q + a x x 2 + a 2 x 2 

2 

3 0 1 3 3 3 

The constants a Q ,a 1 ,a 2 are split into 3 different sets of 
constants corresponding to y^,y 2 ,y^ and then we take the 
average of these to get the final constant. Though one 
might raise serious questions it has been verified that this 
gives ’sufficiently’ correct results. 
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For a first order polynomial this reduces to 



From a set of six -values of errors for six future ranges, 
three such values of a Q and are computed and then the 
average of these values is taken to get the correction fun- 
ction, i.e. 

Error = A Q + Range 

This routine is repeated for calculating the error function 
for range, bearing and elevation. This simplified curve 
fit routine requires 768 bytes of memory storage and has an 
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execution speed, of 5.115 msecs, assuming a time period of 
500 ns. 

A further simplified program to determine the correction 
function is simply averaging the error values. The average 
of the first three values is averaged to ' the next two set of 
three values each and the final result is used as the corre- 
ction function. .Although, this is not as accurate as the 
above curve fit routine this provides a reasonably high impro- 
vement factor in error correction at higher ranges and the 
improvement factor considerably reduces at lower ranges. In 
some cases where the error is marginal, the corrected and 
uncorrected errors are almost identical in this case. However, 
this has the advantage of occupying very little memoiy space 
348 bytes(as compared to 7 r 8 bytes in the previous case) and 
very fast execution time (2.659 msecs as compared to 5.115 msec 
in the previous case). 

In our implementation, the first-order polynomial 
curve fit is preferred to the averaging case because of the 
following reasons. 

.) The corrected error in the first order polynomial curve 
is smaller than the 

fit /corresponding error in the averaging case, at all 
ranges . 

(b) The execution time is not a bottleneck in our case, 
because of the appropriate choice of the sampling 


interval 
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5.4 Memory Requirement 

A major portion of the memory is required for the 
storing of the predicted co-ordinates ,• these being the future 
range, bearing and elevation and the time of flight. The 
factors that affect the size of stack for' each of these are 
the maximum time of flight of shell and the sampling interval. 
For the cases studied, the maximum time of flight encountered 
was 23.7 seconds, however, a slower muzzle velocity or bad 
weather conditions could, hypothetically, increase the time 
of flight to about 40 seconds. Assuming a sampling interval 
of 0.1 second, 400 values of each predicted co-ordinates 
would require storage, jn t succeeding values computed are 
over-written on the previous values as a value once used for 
comparison is not required again. As each value takes 4 
bytes of memory, 1600 bytes arerequired for each predicted 
value. Thus a total of 64OO bytes are required for the sto- 
rage of the predicted co-ordinates. This requirement can be 
halved by reducing the sampling interval to 0.2 seconds. 

Program storage requires 768 bytes for the Main Program, 
768 bytes for the curve fitting routine and 1536 bytes for 
the time of flight routine. 250 bytes are required for 
storing program referenced address esses and 168 bytes for 
storing the constants and other values in fixed memory loca- 
tions. 24 bytes each are required for maintaining the 
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stacks fo- range error, bearing ea.ro r and elevation error. 

This total requirements adds up to 3586 bytes. Therefore, 
the total memory requirement works out to 9986 bytes. Hence 
a memory size of 10K bytes (10240 bytes) is recommended. A 
suggested memory map is given in Figure 5.1 and Figure 5.2 
shows the details of the addresses occupied by the fixed values 
and the errors. 

5.5 Timing Constraints 

The number of time states associated with every 8080 OPU 
instruction used in our program is listed in Table 5.2. 
Instruction execution tine equals number of time periods mul- 
tiplied by the duration of a time period. A time period may 
vary from 480 nanosecs to 2 psec, depending upon the manoiy 
type that is used. For computing the execution time of our 
program - time period -f 500 nanosecs has been assumed. As 
the sampling interval chosen is 100 msec (0.1 sec), the 
number of time periods available for the execution of the 
entire program is 200,000. 

The breakdown of the tine periods required for the 
major subdivisions of the program are listed below : 

Initialisation 479 

Wind components computation 

TIME routine assuming 10 iterations to 
compute time of flight and 5 iterations to 
compute squareroot of a number 


120102 





RBDICT3D FUTURE 
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Stack loading 


504 

Correction routine 


816 

Co-ordinates correction 


516 

t^ storage and clock increment 


277 

clock = t^ check 


. 401 

Error computation 


861 

Stack initialisation 


137 

Curve fitting routine using first 
order polynomial for curve fitting 


10 230 

Retrieval pointer initialisation 


301 

Clock = t„ check 


445 

Mis cell eneous (counter initiation. 

flag check 

237 

etc) 

TOTAL 


135850 


The above figures correspond to maximum time periods for 
each subdivision. This corresponds to a maximum execution 
time of 67.925 millisecs. at start and a maximum of 67 millisecf 
once started (not including initialisation, wind components 
computation and initialisation of time of flight routine). 

This is well within the 100 milliseconds sampling interval, 
i.e. the interval after which the next values of R^, S^, 
are available. 
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Table 5.2 INSTRUCTION EXECUTION TIMES 


MNEMONIC 

NIUBHR OF 'I 

LZI 

10 

PUSH 

11 

POP 

10 

STA 

13 

LDA 

13 

ZCHG 

13 

SPHL 

4 

INX 

5 

MOV r x ,r ? 

5 

MV I r. 

5 

INR 

7 

XRA 

5 

OMP 

4 

JMP 

4 

JZ 

10 

JNZ 

10 

JP 

10 

JM 

10 

DCX 

5 

SH1D 

16 

LHID 

16 

FADD 

100 



contd 


Table 5.2 
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MNEMONIC 


NUMBER OF TIMS 
PERIODS 


FSU3 

FMUIi 

PDF" 

P IF 

oos 

ARSIF 


100 

100 

100 

100 

100 

100 


5.6 System Flowchart 

The complete system flow chart is presented in Figure 5.3. 
The 8080 Assembly Language program corresponding to this is 
listed in Appendix 0. It is assumed that the settings for the 
manual inputs are made prior 'to starting. It is further assumed 
that the wind conditions remain constant. 

The time of flight and predicted co-ordinates are 
obtained from the time of flight routine. This has been 
dealt with in details in Section 2.5. These predicted 
co-ordinates are stored in separate stacks, the pointers for 
which are initialised at the start of the program. A maximum 
of 400 values can be stored. The 401 st. value is 
stored in the 1st location of the particular stack. Thus 
the system will work for maximum initial time of flight of 
40 seconds. The tolerance for checking whether the predicted 
time of flight and clock values are same is restricted to half 




Figure 5.3 
System Flow-Chart 


rhi 
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the sampling interval. This has been found to work very 
satis facte rily for all the ca'-es studied. 

To get early values of the correction coefficients, 
the first six values of error computed are used for the 
curve fit routine and hence the correction coefficients. 
Subsequently, it is not necessary to enter the curve fit 
routine for each subsequent error value as this is redumdanv [ 3 ] 
A prooedure adopted is to enter each 7th value into stack 
and when the stack is full of 6 fresh values, enter curve 
fit routine and obtain fresh correction coefficients. 

Any firing on the target causes an interrupt, and 
the program execution starts again at ’LOOP’. 
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6. OONCLublON AMD SCOPE POP PO&'MBR WORA 

The error correction technique employing first order 
polynomial curve fitting routine has been found to greatly 
enhance the accuracy of the system. Error in range in some 
cases is reduced to as low a value as 0.5 feet while that in 
bearing and elevation is reduced to almost negligible values. 
As the error function ie continuously computed by fitting 
curves repetitively (when the stack of error values is full) , 
fresh correction coefficients are continuously available 
and thus the error values are checked from continuously 
increasing. 

Results also show that a microprocessor can easily 
be adapted to the solution of a real time application. A 
dedicated microprocessor, supported by a small amount of 
specialized interface hardware and a well designed soft- 
ware support is shown to offer a viable alternative to 
special purpose discrete logic system. The low cost of 
microprocessors provide a significant economic advantage 
to systems of this type. 8-bit microprocessor (INTEL 8080), 
with additional hardware, has been found suitable for float- 
ing point arithmetic employing a 24 bits mantissa and 
6 bits exponent. 
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Further work suggested on this topic relates to 

o Simulation of the proposed design on one of the general 
purpose computers and hardware design of the system 
o Implementation of the system on one of the existing 
Fire Control Systems of the Havy . 
o Collecting sufficient data regarding shell trajectories 
through a number of experiments and trials and compiling 
tables which relate gun elevation and range to the time 
of flight. With these tables available studying the space 
■ . . -time trade off in calculating the time of flight by the 
iterative method used in our design and that by using 
Table look-up and interpolation, 
o With the present availability of a wide variety of 16 bits 
microprocessors, investigating the implementation of this 
design using a 16 bit microprocessor and eliminating 
the need of support hardware for floating point arithmetic. 
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APP3EDTZ A 

ffeperation Corrections 

Jxl tiie text it was ass umed that the gun and the radar 
antenna are theoretically at the sane point and plane. How- 
ever, this is not so practically. To account for this 
separation of the gun (G-) and director (D) the undermen- 

"I 

tioned corrections are applied. 

(a) DIP Gorrection : The correction to elevation necessary 
because of the difference in height above deck-plane of 
director and gun. It can be defined as the angle between 
the director line of sight and the gun line of sight due to 
the vertical displacement (h) between the gun and the dire- 
ctor. 



Sin Dip = SG/% = h. Cos (S f )/H f 
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, b) Co nver g ence Girrection : In h.A. fire this is required, 

in both the horizontal and vertical planes. This is due to 
the horizontal displacement (d) between the director and the 
gun, Diagracuaati cally these are explained below ; 




Sin 0 = SG/E f = d Cos B f Sin S f /B f 



PREDICTION ERROR PROGRAM 


PROGRAM IT. COMPUTE ERRORS IN PREDICTION 


DIMENSION TOF(20Q),pRF«2QQ>,PBF<2QO)#PSF<2QO>#PRFC<-20Q) ,PBFCC2O0), 

1PSFC (200 ) 

COMMON RP , SP, BP, RRA TE,SRATE, BRA TE,DT,TF, VI, SF,8F,RRF,DTR,L,M 

common g, jstart»wr,wl 

COMMON AA(3),BB<3),STAC(6,3},RANGE<6>,RFA<200)»SFA(2Q0)#BFA<200) 

INTEGER P 
DTR * 22 ,/1260, 

VI * 2350, 

DT » 0,2 
Q s 32, 

HRIT6<4,1001 . . . .... 

100 FORMAT (//#—-**- — — — — — — 

2//#DIMENSlON FOR PARAMETERS AR£#,/*RAMGE IN FEET*,9Xi*RANGE RAT 
3E IN FEET PER SECOND#, /*ELEVATIQN IN DEGREES*, 2X,*ELEVATI ON RATE I 
4N DEGREES PER SECOND*, /*BEAR ING IN DEGREES#, 4X , *BEAR ING RATE IN DE 
5GREES PER SECOND# ) 


WRITE <6, 101) . . \ . 

101 FORMAT </*TIME IN SECONDS*, /*WIND VELOCITY(W) IN FEET PER SECOND*, 
1/#WHD INCLINATION IN BEGREE$(WI )*,/*RANGg ERROR IN FEET*, /*ELEVAT 
3 I ON and BEARING errors IN minutes#,//*— — ............ 


Do HU I * 1,14 


>--#///) 


RP, SP, BP, ARE THE PRESENT RANGE , ELEVATION , AND BEARING RESPECTIVE^ 
RRATE,SRATE,8RATE,ARE THE COORESPONDING RATES \ 

DT IS THE SAMPLING INTERVAL 


WIND VELOCITY IS 0 TO 100 FEET PER SECOND 

WIND INCLINATION IS ASSUMED 10 DEGREES TO 70 DEGREES 



RE AD ( 5 , 10 ) RP,SP, BP, PRATE, SRATEiBRAJE 
10 FORMAT(F10,2,F7,2,F8,2,F8.2,F7,2,F7.2) 

RP X R P # 3 '» ■ ' . ■ . - 

NUMW * pUMW + 1 

W a IOC , *RNDY1(DUM1) 

WID * 10, + 60, * RNDY2(DUM2) 

WRITER, HO) NUMW , ■; 

110 FORMAT (/X# INPUT DATA , *,I2,*,#) . 

WRITE(6, 103)U, WID, RRATE,SRAt8, IRATE, RP, SP, 8p . , . . . 

103 FORMAT (#CONSTANT PARAMETERS ARE*,3X#* W x #,F12,5,3X,*)<iI * *#F12,,!5 
t,3X,*RRATE * *»F8,3,3X,*SRATE = #,F6,3,3X,*BRATE »*,F7,3,/*INITIAL 



2 PARAMETERS ARE* , 3X» *RP a *,Fl2,3, 3X,*Sp = *,F7,3,8 x,* BP = », 

3F7 i o ) 

hi = KID « D T R 
nR 3 w # COS ( W! ) 
k : L = W * SIN < NI ) 

I START s i 
P s G 

JSTART S o 
ITRIG = 0 
N 3 u 
CLOCK = C 
H s % 

L = G 

I COUNT s 0 
LCGUNT « 0 
K s 0 
Hi a 1 

M2 * 1 

FIREAT * ^000, 

RELAX » 2* ft RP 

IF ( RP ,LT » 20001 1 ) RELAX = 3, « RP 

SP * SP * DTR ' 

BP ® BP * DTR 

BRATE * BRATE#DTR 
SPATE a SRATE*DTR , , . 


L IS A FLAG TO DETERMINE IF THE CLOCK VALUE CORRESPONDS TO A VALUE OF 

COMPUTED TIME OF FLIGHT 


20 IF C L ..EG* 1 ) GO TO 40 
N = N + 1 
CALL TIME 

IF(JSTART ,EG, 1 ) GO TO 11 

JSTART * i 
WRITE (6,104) TF 

104 F0RMAT(1X,»INITIAL TIME OF FLIGHT s *,FlO,4) 

11 IF ( RP ,LE, FIREAT ) ITRIG a 1 
IF C ITRIG , EG, 1 ) GO TO 1110 
IF ( RRF if EQ., 0, ) GO TO 1111 
IF ( RP ,GE , RELAX ) GO TO 1110 
PRF ( N ) a RRF 
PSF(K) s SF 
PBF(N) a BF 

RFA(N) s RP + RRATE « TF 

DV a ARSIIV ( RP / RFA(N) * SRAT& * TF ) 

SFA(N) ■ SP + DV 

DA = ARSIN.(Rp*C0S(SP)*BRATE*TF/(RFA<N)*C0S(SFA(N) >)) 
BFA(N) a BP ♦ DA ' ' ' . , 

IF (BFA (N ) ,GT, 6,28 ) BFA(N) a, BFA (N ) • 6,28 
I F ( BF A ( N ^ ,LT, 0, ) 8F A ( N ) a 6,28 + BF A ( N > 

IF (K ,EQ, 1 ) GO TO 22 
PRFC(N) a PRF(N) 
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P&FCCiv 5 a PSF(N) 

P 8 F C C |V ) a P3F(N) 

GO TC 21 

22 CALL C0R6CT ( PRF ( N ) , DRX ,DSX , DBX ) 
PRFC(fv) a PRF !N ) - DRX 
PSFC!N) « PSF(N) - QSX 
P8FC (N ) a PBF(N) - DBX 
21 TOF < N ) » CLOCK + TF 
CLOCK1 * CLOCK * DT 
30 0 IFF a ABSCCL0CK1 . • TOF(M) ) 

IF ( DIFF ,LT ♦ (DT/2’,) ) GO TO 31 

CLOCK a CLOCK1 

GO TO 32 


IF L * 1 THEN THE VALUES OBTAINED FROM INPUTCRADAR) ARE COMPARED WITH 
THE COMPUTED VALUES AT THAT INSTANCE TO DBTAIN THE RELEVANT ERRORS, 


V * *#» , 

IF < ISTART , EQ , 0 ) GO TO 20 
ISTART * 0 

32 CALL COMPUT(CRF,CSF, CBF) 


COMPUT RETURNS THE FUTURE POSITION COORDINATES AFTER SAMPLING INTERVAL 

TIME ADVANCED BY SAMPLING INTERVAL DT 

THE PRESENT COORDINATES NOW BECOME THE COORDINATES AFTER SAMPLING TIME 


RP = CRF 
BP a CBF 
SP = CSF 

IF ( N ,ES, 200 ) N a 0 

GO TO 20 

40 CALL ERRORS! RE»EE,8E, PRF <M),PSF<M),PBF CM)) 

DR a RE 
DE a EE 
DB a .BE 

DEM * DE / DTR * 60, 

DBM a D8 / DTR * 60, 

IF C K ,E0. 0 ? .GO .TO. 41 . ........... 


WHEN K a l CURFIT ROUTINE HAS BEEN ENTERED HENCE CORRECTED PARAMETER 


CALL ERRORS! RE,EE # BE,PRFC(M),PSFC<M)#PBFC<M)> 

DRC a RE 
DEC a EE 
DBC a be 

DECM « DEC / DTR * 60, 

08 CM « DBC f DTR * 60. 

IF ( DR ,EQ» DRC ) GO TO 41 



LCCUM » LCOU^T + 1 

iraccuM .eg, lo > lcount * 1 

IHLCOUNT »N£» 1 ) GO TO 41 
WRlTE<6»10S)CLOCKl»DR,DRC # DEM»DECH»DBH*DeCH 
105 FORKAT</«AFTER #,F10,4»* SECONDS OF DETECTION *, 

1/# UNCORRECTED RANGE ERROR * *,F11, 5,lOX,*CORRECTED = *,Ftl',5, 

2/* UNCCRRECTED ELEVATION ERROR * * ,FU,§flOX, "CORRECTED * *,FU,5, 

3/# UNCORRECTED BEARING ERROR * *,F11 » 5, 1QX, "CORRECTED s *,F11»5) 

I COUNT * I COUNT + 1 
IF! I COUNT ,EQ, 6 ) GO TO 1110 
41 L * a 

IF ( HI ,EG, 6 > GO TO 42 
IF { P .EG, 0 ) GO TO 43 

GO TO 44 . , , . , 


P M FOR H ,LE, 9 AND THEN IS ALWAYS 1 


42 P ■ 1 
Ml. • 1 . 


Ml IS A COUNTER THAT DECIDES ENTRY INTO DTACK 


43 RANGE ( M2 ) * PRF(M) 
STAC(M2,1) = DR 
STAC ( M2 * 2 ) = DE 
STACIM2 »3 ) a DB 
M2 a M2 + 1 


M2 KEEPS TRACK OF NUMBER OF ENTRIES IN STACK 


IF { M2 ’,EQ» 7 ) GO TO 45 
44 M a M* 1 

IF ( M »EG, 201 ) M * 1 
Ml * Ml + 1 
GO TO 30 


45 M2 * 1 



STACK POINTER INITIALISED 


CALL CURFIT 
K a 1 
GO TO 44 

110 WRITEUf 1C6)CLQCK1*TF , 

106 F0RMAT</*FIREQ AT TARGET AFTER *»F10,4»* SECONDS OF DETECTION. 
1TIME OF FLIGHT a *»F10,4) 

111 CONTINUE 
STOP 
END 



T I We ROUTINE 


S U O H C 1 T I i\ £ TIKE 

CO?MN oPs?AR?"wK E ' SBAT6 ' eRATE ’ ET ’ TF ' V1 ' SF - BF ’ RR,r ' D7R '‘-.« 
RSrL C L«!NMR B<3, ’ ST4C<6 ’ 3, ’ RANGE!6, ' RFA<2001 ' SF -'< 2 00>*BrA(200) 


IF ( JSTART f EG, i ) go TO 20 

v * VI 

OLDR s RP 
OLDT * RP / V 
SF a SP 

20 A * G # S IN (SF ) 

JSf-SS * LT # . 45000 * ) \ 

IF ( OLDR ,GT, 60000,) 

B a -2, * V 
c a 2, # OLDR 

XI 9 0 # 

X2 a 0 i 

CALL ROpTS, <A,B,p,Xl,X2) 
IF((X1 ^QT,, 0,) ,♦ AND., (X 
IF((X1 ,GT , 0,J , AND, (X 
GO TO 44-4 

30 NEWT * XI 
GO TO 31 

32 NEWT a X2 

31 TERM1 « RRATE • NEWT 


V = VI 


10, * Vi 


,6 and,. 

.AND, 


0 ,)) 

0,)) 


GO TO 30 
GO TO 32 


0,05? GO TO 40 


- OLDR) ;GT, 1,00) Go TO 40 


NEWR a RP + TERM1 

IF (ABS (NEWT - OLDT) *,GT, 0,05? GO 
IF ( A8S ( NEWR - OLDR) , GT, 1,00) Go 

GO TO 50 " 

40 N s N+l 

IF ( N .EG, 15) GO TO 555 

OLDT = NEWT 
OLDR s NEWR 

EXP1 = RP/NE.WR * SRATE *nEWT 
IF (EXPl ,GT, 1 , ) ExPl = 0,98 
IF (EXpl ,LT , (-l # ) ? EXP1 a (*,98 ) 

DV * ARSIN(EXPI) 

DVD a 0V/DTR 
SF a SP + DV 

IF ((NEWR jLT.t 30000’, ) AND, (SF ,LT 
IF ( (NEWR , LT , 30000,) .AND , (SF .GT, 
IF (SF ,LJ, 0, ) GO TO 44 
IF (SF ,GT, 1,571427) SF a 1,54 

41 SFD a SF/DTR 
GO TO 20 

42 SF a *,174603 

GO TO 41 

43 SF a l', 3619034 
GO TO '41 

44 IF (NEWR ’,IT, 60000, ) SF a 0.0174602 
SF a SF * DV 

GO TO 41 

444 SP * SP ♦ 5,#DTR 


.AND, (SF 
.AND, (SF 


» 174603)) GO T(? 42 
1,396824) ) GO TO 43 


) SF a 0,0174603 



555 WRITE<6»556) 

556 Format (img»iox#«lqst target #> 

RRF a 0» " ' 

return 

50 EXP1 * (RP#COS(SP)*BRATE«NEWT )/(NEWR *CGS ( SF ) I 
IF (EXP1 ,, GT, 1, > EXP1 » 0,98 
IF ( EXP1 , LT , (-1, ) ) EX PI a <-,98 > 

DA « ARSIN(EXPI) ' ' 

BF a Bp * DA 
BFD * BF/DTR 

IF (BFD '.-UT « 0, }, BFD = 360, + BFD 
IF ( BFD ,GT, 360, ) BFD = BFD - 360, 

BF * BFD # DTR 
TF a NEWT 

RRF s NEWR . . 


INTRODUCING CORRECTIONS 


TANEUE » G * TF # COS(SF) / ( 2, * V) 

PHI a SF ♦ TANEUE 
DRICON a 0,04 

drift a dricon « tanele 

BF a BF - DRIFT 

DEtV = O', 10 
DEIC = 0,10 , 


CORRECTIONS DUE TO CHANGES IN NUZZLE VELOCITY AND BALLISTICS 


IF (JSTART ,EQ, 1 ) GO TO 203 

DELTM = RNDYl(DUHl) # DELV + RNDY2 ( DUM2 ) * DELC 

DELDVH a (RNDY2(DUM2)*PELV + RNDY1 (DUHl) ft DEUC ) ft . 0 ,1 ....... . 


effects due to wind 

CHANGE IN RANGE,ELEVATION AND BEARING DUE TO WIND 


203 DELRW « WR ft TF ft COS ( PHI ) 

DELPHI a - ( WR * TF ft SJNt PHI ) ) / RRF 

DELAW a ( WL * TF ) / < RRF ft COS < SF ) ) 


THEREFORE FINAL FUTURE RANGE ELEVATION AND BEARING ARE 


204 RRF a RRF ♦ DELRW 

GUNELE a PHI * DELPHI 
8F a BF + DELAW 



no on 


IF ( BF ;.GT, 6,2857) BF » BF • 6,2857 
IF<BF t LT, 0, > BF s BF * 6,2857 


THE RANGE CORRECTION IS EFFECTED BY CORRECTION TO TIME OF FLIGHT 

AND ELEVATION AS FOLLOWS 


DELTW * HR « TF » C0?(SF) / V 

DELOVW «* «<(NR # COS(SF) ) / < V * SIN(PHI) * COS(PHI) ).), 
J THEREFORE FINAL TF AND SF ARE 


205 TF * TF* ♦ DELTW + DELTM 
SF » SF + DELDVW + DELDVM 

RETURN 

END 

SUBROyTl^e ROOTS (A#8,C,X1*X2) 

IF (A , EG , P,>. RETURN 

ir<e#B , lt * 4, *a*o, Return 

Xi ¥ OB fSQRT(B*B ? (4, t *A*C) ))/(2,#A) 

X2 * OB *$GRTCB*BM4,*A*C)))/<2»#A> 

RETURN 

END 


AIRCRAFT SIMULATION ROUTINE 


SUBROUTINE CQMPUT(CRF,C$F,C8F) 

COMMON RP ,SP # 8P,RRATE # SRATE,8RATE,DT»TF,V1,SF,BF»RRF,DTR,L#M 

COMMON G , USTART , WR* WL , , , ... 

COMMON AA(3),8B(3),STAC(6,3),RANGE(6)#RFA(2OO)*SFA(2G0>#BFA<2OQ) 

TERM1 * RRATE * DT , 

TERM2 a O', 5 # RP*DT *DT 

TERM3=SRATE#SRATE4BRATE*BRATE#COS(SP)#COS(SP) 

TERM3 » 0, 

TERM2 » TERM2 « TERM3 
CRF * RP * TERMl * TERM2 
DV ■ ARSIN <RP/CRF*SRATE«DT) 

CSF a SP ♦ DV 

DA s ARSIN ( RP*COS<SP)*BRATE»DT/(CRF*COS<eSF)) > 

CBF a BP + DA 
CBFD b CBf/DJR . 

IF (CBFD ,LT, 0, ) CBFD. 3 360, * CBFD . 

IF (CBFD ,GT, 360, ) CBFD a CBFD * 360, 

CBF s CBFD 3 DTR 

RETURN 

END 


10 0 


ERRORS ROyTlNE 


SUBROUTINE ERRORS (RE, EE, BE, PRF,PSF.PBF) 

COMMON RP ,Sp,BP,RRATE,SRAT£,BRATE,DT,TF,Vl»SF,8F , RRF , DTR#k*M 

COMMON Q,gSTART,WR,WL 

COMMON AA(3),8B(3),STAC(6,3),RANGE<6),RFA(2QQ)-,SFA(20Q),BFA<2QQ) 
RE * PRF * RFA(M) 

EE » P$F - SFA(M) 

BE * PBf * BFA<M> 

IF (BE , QT ,, 3,14 BE a 6,28 • BE 

IF ( BE ,LT, ( -3,14 ) ) BE = 6,28 * BE 

RETURN 

end 

ibftc subs 


CURVE FITTING ROUTINE 


SUBROUTINE CURFIT 

COMMON RP ,SP,BP,RRaTE,SRaTE,BRaTE,DT,TF,V1,SF,BF,RRF,DTR,UM 
COMMON G,USTART,WR,WU 

COMMON AA(3),BB(3),STAC(6,3) ,RANGE(6),RFA(200),SFA(200),BFA(200) 
DIMENSION AO ( 3 ) » Al < 3 ) 

L * 1 
10 I « 1 

J s 1 

DO 1 K a 1,3 
XI = RANGE! I) 

X2 = RANGE ( 1+1} 

Yi » STAC ( I it ) 

Y2 * $TAC(I + 1,U 
D a X2 » Xl 

. AO < J ) * (X2 * Yl - XI * Y2 ) / D 

Al(U) » ( -Yl + Y2 ) / D 

J = U ♦ 1 
I a I + 2 

1 CONTINUE 

AA(L) * (A0(1) + A0(2) f AO <3 > ) / 3, 

8B(U a (Al(l» ♦ Al( 2) + Al(3* ) / 3 f 

L . ® . 4* % 

IF ( L . LE, 3 ) GO TO 10 

RETURN 

END 
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APPENDIX C 

I MACROS USED IN THE ASSEMBLY LANGUAGE PROGRAM OF THE SYSTEM 

««««««## ««««*« 

f macro name is in the label field 


TRANS 

MACRO 

L0C1,L0C2 

.MACRO TO TRANSFER 4 BYTES 
{ FROM MEMORY LOCATION. LOCI 
ITO MEMORY LOCATION L0C2 


LHLD 

LOCI 



SHLD 

L0C2 

• . ' ■' ‘ ”* 


LHLD 

LOCI ♦ 2 



SHLD 

L0C2 * 2 



ENDM 



INSRT 

MACRO 

SPX 

J MACRO TO INSERT A 4 BYTE 
} WORD INTO STACK 

ISPX indicates current val 
{ OF stack pointer 


MOV 

B # H 



MOV 

C #L 

' ' ■ ; ' ' ■ ' : 


LHLD 

SPX 



SPHL 


■ : , ' ■ . ■ . ■ ■■ '■ : ! 

■■ ■ . '■ i 


PUSH 

D 

. ' ■ ' • ( | 


PUSH 

B 

{STACK POINTER ADJUSTED 


DCX 

H 


DCX 

H 

i TO INSERT NEXT VALUE ! 


DCX 

H 



DCX 

H 



SHLD 

SPX 

■ . " J 


ENDM 



RETRY 

MACRO 

IPX 

j MACRO TO RETRIEVE A 4 




{BYTE WORP FROM STACK 
# IPX IS POINTER TO BYTE 




{TO BE REMOVED 


LHLD 

IPX 



SPHL 




POP 

H 



POP 

D 



ENDM 



DETER 

MACRO 

ADDR1, ADDR2 

{MACRO TO DETERMINE THE ! 
{DETERMINANt TO FIND A 



CCEFF 


LXI SP# ADDR1 

POP H 

POP D 

SHLD XI 

XCHG 

SHLD XI t 2 

LXI SP# ADDR2 

POP H 

POP D 

SHIP X2 

XCHQ 

SHLD X2 ♦ 2 

XCHG 

FSUB XliD 

SHLD DET 

XCHG 

SHLD DET ♦ 2 
ENDM 


MACRO ADDR1,ADDR2,STQ1#ST02 


LX I 

SP # ADDRl 

POP 

H 

POP 

D 

FMUL 

X 2 # D 

SHLD 

TEMP 

XCHG 

SHLD 

TEMP + 2 

LXI 

SP* ADDR2 

POP 

H 

POP . 

D 

FMUL 

XI, D 

FSU8 

D*TEMP 

FDIV 

det,d 

SHLD 

STOl 

XCHG 


SHLD 

STOl + 2 

LXI 

SP# ADDRi 

POP 

H 

POP 

D 

FSUB 

D, ADDR2 

FDIV 

i i ' r*s 

DET,D 

■ mm mm 


{ solution for obtaining 

I COEFFICIENTS OF FIRST . 

I ORDER POLYNOMIAL 
i ABDRli ADDR2 REFER TO TWO 
{ RANGE VALUES 


| S' 


{ MACRO TO COMPUTE THE TWO 
{COEFFICIENTS OF A FIRST 
{ORDER PILYNOMIAL i 

1X1, X2 ARE TWO RANGE VALUE ! 
{STACK HAS ERROR VALUES . 
{COEFF', STORED IN ST01,ST02 
{ADDR1,ADDR2 ARE POINTERS 
{TO ERROR VALUES 




k 

I 

\ 


I 

■■ 


XCHG 

SHLD 

ENDM 


ST02 * 2 


AVRGE MACRO LQC1,L0C2,L0C3,LQC4 {MACRO TO DETERMINE THE 

I final COEFF, BY averaging, 
13 SETS OF COMPUTED COEFF, 
5 LOCI | L0C2', L0C3 CONTAIN 
{THE VALUES TO BE AVERAGED 
JL0C4 STORES FINAL COEFF, 

LOCI 4 2 

HOC1 
LQC2,D 
L0C3d D 

#D ID IV IDE BY 3 

L0C4 

L0C4 + 2 


decrt 

MACRO 

I P T 

{MACRO 

TO 

DECREMENT 

THE 




{VALUE 

OF 

A POINTER 

BY 4 


LHLD 

I FT 






DCX 

H 






DCX 

H 






DCX 

H 






DCX 

H 






SHLD 

IPT 






ENDM 






RTVAL 

MACRO 

ARGl, ARG2 

{MACRO 

TO 

STORE THE 

RIGHt 


{VALUES OF PREDICTED CO<*RD 
{ARG1 CORRESPONDS TO. THE. 
{COORDINATE TO BE CORRECTED 
{AND ARG2 TO THE CORRECTlO 
{FACTOR, 

LHLD ARG1 *2 

XCHG 

LHLD ARG1 • ( 

FSUB ARG2# D 

SHLD ARG1 

XCHG 
SHLD 
ENDM 


LHLD 

XCHG 

LHLD 

FADD 

FADD 

FDIV 

SHLD 

XCHG 

ShLD 

ENDM 


ARG1 + 2 


» THE CORRECTION COEFF(STO) 
CGEF1, C0EF2 ARE THE COEFFS 
?OF THE EQUATION 
{ERRORsCOEF1+RANSE»COEF2 

LHLD NEWR ♦ 2 
XCHS 

UHLD NEWR 

FMUU DiC0EF2 

FADD COEF1 

SHUD STO 

XCHG 

SHUD STO + 2 

ENDM 


f ######## ■»#*#•##* <Mt -iHUHtfl* {HHHHHHHtttftlHHHl 

JMAIN PROGRAM IN ASSEMBLY LANGUAGE OF 6080 


}««*•# #* #* #* #* *# #* ** ** ## *# ** *» *« *«■ »#■ *# ** #* *# «« #* ## «« *« «# ## «««« 


J PROGRAM ASSUMES RP # SP,BP,R,S,B, ARE CONTINUOUSLY AVAILABLE 
I IN MEMORY LOCATIONS SPECIFIED FOR THESE, 

^SETTINGS FOR H*HI iV1*D£LVjCELC ARE TRANSFERRED to MEMORY 
f LOCATIONS RESERVED FOR THESE 


\ 


j *#*###** *##**#**#*###*#***####«*****#**#***##*###*****»#****# 

J INITIALISATION 

i *■&#**#*«■»■»*•###•»■»«•**#•*■&###«••#*#*###*#■»«**###*#**»*#**«***«■***** 


STARTS XRA 

A 

t VALUES INITIALISED 

STA 

N 

. i TO ZERO 

STA 

USTAR 


, STA 

P 


STA 

K 


STA 

L 


MOV 

H » A 


MOV 

L » A 


SHLD 

CLOCK 


SHLD 

CLOCK + 2 


MVI 

A# 1 

i VALUES INITIALISED 

STA 

M 

1 TO ONE 

STA 

Ml 


STA 

M2 


LXI 

Hi 2761H 

t POINTERS. INITIALISED 

SHLD 

SP1 

SSTORAGE POINTER FOR PRF 

LXI 

Hi 2121H 


SHLD 

SP2 

• J STORAGE POINTER FOR PSF 

UI 

H.1AE1H 


SHLD 

SP3 

^STORAGE POINTER FOR PBF 

LXI 

H»14A1H 


SHLD 

SP4 

SSTORAGE POINTER FOR TOF 

LXI 

H,275DH 


SHLD 

IPT1 

|RET REIVE POINTER FOR PRF 

UI 

H,2UDH 


SHLD 

IPT2 

IRETREiVE POINTER FOR PSF 

LXl 

h.iaddh 


SHLD 

IPT3 

SRETReIvE POINTER FOR PBF 

LXI 

H,149DH 




SHLD 

IPT4 

{RETREIVE POINTER FOR TOF 

L X I 

H,0E61H 

SHLD 

STAC1 

{POINTER FOR DR 

LXI 

' H| 0E49H 

SHLD 

STAC2 

{POINTER FOR DE 

LX I 

H,0E31H ■ 


SHLD 

STAC3 

{POINTER FOR DB 

LXI 

H,0E19H 

SHLD 

STAC4 • 

{POINTER FOR RANGE 1 

COMPUTE WIND COMPONENTS 


EQUIVALENT WIND VELOCITY(W) IN LOCATION 0D1DH 

WIND 

INCLINATION! WI ) IN LOCATION 

QD19H 

LHLD 

XCH6 

0D18H 


LHLD 

cos 

0D19H 


FMUL 

D,0D1DH 


SHLD 

NR 

{ WR a K * COS c HI > I 

XCHG 


SHLD 

WR ♦ 2 


LHLD 

XCHG 

OD1BH 


LHLD 

0D19H 


SIN 



FMUL 

DiODIDH 


SHLD 

WL 

{ WL * W * SIN(WI) 

XCHG 



SHLD 

WL + 2 



J CHECK IF ERROR FLAG IS SET 

LOOP! LDA L • ' 

JNZ ERORS j FLAG SET IF NOT ZERO 

LDA N 

I NR A 

MV I '8,401 

CHP B 

vlNZ CQNTN ? CONTINUE if N ',LE, 400 

I N » GT » 40O # THEREFORE INITIALISE ALLSTORE POINTERS AND N 

LX I H,2?61H 

SHLD SP1 

LX I H.2121H 

SHLD SP2 

LX I H,1AE1H 

SHLD SP3 

LX I H,14A1H 

SHLD SP4 

MV I A , 1 | 



: time returns the time or flight 

and predicted coordinates 


MV I 

A#1 


* 1 


S.TA 

USTAR 



TIMER: 

LHLD 

NEMR + 2 


ICHECK NEWR 


MOV 

A.#H 


i IF NEWR * 0 TARGET IS LOS 


JNZ 

KCHEK 




XCHG 





LHLD 

NEWR 




MOV 

A i H 




JNZ 

KCHEK 




UMP 

START 



KCHEK: 

LHLD 

NEWR + 2 


, pur PREDICTED COORBDS, | 


XCHG 



* IN THEIR RESPECTIVE ] 


LHLD 

NEWR 


; STACKS 


INSRT 

SP1 


' 


LHLD 

SF + 2 


' ■ " : 


XCHG 





LHLD 

SF 




INSRT 

SP2 




LHLD 

BF + 2 




XCHG 





LHLD 

BF 




INSRT 

SP3 




LDA 

K 




JZ 

PRDCT 




JMP 

CORCT 


# STORE CORRECTED PREDICTED 





jvalues 


l CORCT RETURNS THE 

CORRECTION FACTOR 

C0RCR5 

RTVAL 

NEWR,DRX 


l 


rtval 

SFiDSX 




rtval 

8F.DBX 



PRDCT : 

LHLD 

CLOCK + 2 


J increment by value OF Ci 


XCHG 





LHLD 

CLOCK 




FADD 

D,NEWT 




INSRT 

SPA 




LHLD 

CLOCK 




FADD 

D,0D99H 




SHLD 

CLCK1 




XCHG 



1 CLCKZ 5 CLOCK ♦ DT 


SHLD 

CLCK1 + 2 




1 CLOCK1 

AND TOF(M) 

compared 


TLOOP: 

RETRV 

IPT4 




F SUB 

d,clcki 




MOV 

AiD 




UP 

ABSOL 




TRANS 

CLCK1, CLOCK 


JMP 

ILOOP 

1 

LFLAGJ H VI 

A#1 


STA 

L . . . 

1 

JMP 

ILOOP 


{ ##**###*#*#**&**#*#*#*#*##«##**###«**####*######*.##*«*#*#***# 

{ ERROR 

calculation routine 

■ | 

' 

- ■ ^ ; 


ERORSS XRA 

A 

■ 

STA 

L 

. 

L DA 

Ml 

■ ' 

MV I 

B,6 

LOAD FIRST 6 VALUES IN STA j 

CMP 

B 

SAND THEN EVERY SIXTH VAL, 

JZ 

resto 


LDA 

P 


JZ 

PILE 


JMP 

I NORM 


RESTOt MV I 

A, 1 


STA 

P 


STA 

Ml 


PILES RETRV 

IPT1 

{RANGE FOR WHICH ERROR IS 

INSRT 

STAC4 

; COMPUTED IS STORED IN STA 

MOV 

Hi8 


MOV 

L# C 


psub 

OD01H »D 

{ERRORS COMPUTED AND STORE 

INSRT 

STAC1 

1 in stack : 

RETRV 

I PT2 


FSU8 

0D05H » D 


INSRT 

STAC2 


RETRV 

IPT3 


FSUB 

OD09HiQ 


INSRT 

STAC3 


LDA 

A, M2 


INR 

A 


MV I 

8,7 


CMP 

B 


JZ 

STAIN 

{STACK FULL’, INITIALISE 

STA 

M2 


JMP 

INCRM 


STAIN: MVI 

A , 1 


STA 

M2 


STA 

K 

{CORRECTION FLAG IS, SIT 

LXI 

H.QE61H 

{INITIALISE STACK POINTERS 

shld 

STAC1 


LXI 

H, 0E49H 


SHLD 

STAC2 



I CURVE FITTING ROUTINE 


t ««««»« #«*«#« ##*##** «######*# 4*# *##*######*# a##***### #*«##*##* 

DETER 0E15H* OEijH iGET THE DETERMINANT 

COEFF OE5DH,OE59H,AORl,AlRl 

COEFF 0E45H»0E4iH,A0El,AlEl 

COEFF 0E2DH,0E29H,AQB1,A1B1 

DETER OEODH, 0EQ9H 

COEFF QE55H,QE51H*A0R2,A1R2 

COEFF 0E3DM,QE39H,AQE2,A1E2 

COEFF OE25W,0E21H,AOB2,A1B2 

DETER OE05H * 0E01H 

COEFF OE4DH,OE4$MOR3#AlR3 

COEFF 0E35H # GE3lM,AQE3,AlE3 

COEFF QE1DH,QE19H»AGB3,A1B3 

AVRGE AORi,AOR2,AOR3,AAR 

AVRGE AOE1,AQE2|AOE3,AAE 

AVRGE A0B1«APB2 # A0B3 # AAB 

AVRGE A1R1»A1R2# A1R3,BBR 

AVRGE A1E1»A1E2,A1E3 # BBE 

AVRGE A1B1,A1B2»A1B3,BBB 

JEND OF CURVE FIT ROUT INE, INCREMENT RETRIEVE POINTER M. 

{DECREMENT RETREIVAL STACK POINTERS BY 4 


RMt IDA 

M 

I NR 

A 

MV I 

B.4Q1 

CMP 

B 

JZ ' 

INITM 

STA 

M 

IDA 

Ml 

INR 

A 

STA 

Ml 

DECRT 

IPT1 

DECRT 

IPT2 

DECRT 

IPT3 

DECRT 

IPT4 

UMP 

TLOOP 

Th * MV I 

A»1 

STA 

M 

LXI 

H, 275DH 

SHLD 

I PT1 

LXI 

H# 211DH 

SHLD 

IPT2 

LXI 

H ilADDH 

SHUD 

IPT3 

LXI 

H»149DH 

SHLD 

j P T4 

JMP 

TLOOP 


{#«#####»#*##*###*#**####**#»*#*#«*******«*#»****#*»********■** 


u CORRECT ROUTINE, THIS DETERMINES THE CORRECTION FACTORS, 




CORCTI FUNC 
FUNC 
FUNC 
. UMP 


AAR#B8R f DRX 
AAEfBBEfDSX 
AAB »8B8*fiBX 
CORCR 


}«•«*««#«»«*«#««##«#«*«««««* «««««««« ««««## ««#««»«« | 
I TIME ROUTINE 1 

|«#«#«# ««««#« «««#«««*##«««««««#«#»# ««#««« «« #««««« ««««#« #««««»« | 


I THIS COMPUTES THE TIME OF FLIGHT AND THE FUTURE COORDINATES 


TIME? 


LOOP: 


INCRVI 


MV I 

A, 15 

STA 

NN 

LDA 

JSTAR 

UNZ 

LOOP 

TRANS 

QD25M,V 

TRANS 

ODOlHfOLDR 

TRANS 

0DQ5H #SF 

LHLD 

0D27H 

XCHG 


LHLD 

0D25H 

FD IV 

D, ODOlH 

SHLD 

oldt 

XCHG 


shld 

OLDT * 2 

LHLD 

SF * 2 

XCHG 


LHLD 

SF 

SIN 


FMUL 

D* 0D55H 

shld 

AA 

XCHG 


SHLD 

AA ♦ 2 

iCHECK 

IF RANGE ,GT 

(TO FIX 

MUZZLE VELOI 

LHLD 

OLDR + 2 

XCHG 


LHLD 

oldr 

FSUB 

D#0D91H 

UM 

incrv 

TRANS 

0D25H » V 

UMP 

calcb 

LHLD 

OLDR ♦ 2 


6 0.0 00 FEET OR 


i V « VI 


l AA * G * SIN(SF) 
LT i 45000 FEET 


IV = VI IF ',LT* 45000 FEET 


XCHG 



LHLD 

QLDR 


FSU8 

D» 0D95H 

{ 60000 - OLDR 

UP 

CALC8 


LHLD 

0D27H 


XCHG 



LHLD 

0D25H 


FMUl 

D# 0D89H 

J VI * 10 

• SHLD 

V 

XCHG 



SHU 

V + 2 

? V 5 1Q*V1 IF OLDR »GT« 6| 

CALCBl LHID 

V * 2 

XCHG 



LHLP 

V 


F HUL 

D,0D79H 


SHLD 

BB 

5 68 = *2 * V 

XCHG 



shld 

88 + 2 


LHU 

OLDR + 2 

1 

XCHG 


E 

LHLD 

OLDR 

' ' '' ' ' ^ ^ ' , | 

FHUL 

D# 0D75H 

: ' , ' ' ' ' | 

SHLD 

CC 


XCHG 


' ' ' ' ^ ' | 

SHLD 

CC + 2 

{ CC a 2 # OLDR ; 

UMP 

ROOTS 


? ************************************************************* 

# ROOTS RETURNS VALUE OF 

TIME OF FLIGHT IN NEWT 

# ************************************************************* j 

ROOTRJ LHLD 

NEWT + 2 

' ■ ' ■ ' ■ '' . ' 

' '' , j 

XCHG 

LHLD 

NEWT 


FSUB 

D * 0D6DH 

{ 0 - NEWT 

MOV 

A i D 

, . " <h . • J *% . , ' f-v 

UP 

L00P2 

I IF NEWT « 0 » TRY AGAIN 

I ALTERING SF 

LHLD 

NEWT + 2 


XCHG 



LHLD 

NEWT 


fmul 

D , ODQDH 


FADD 

D» CDQ1H 


SHLD 

NEWR 

JNEWR » RP ♦ R # TF 

XCHG 

SHLD 

NEWR + 2 


{COMPARE 

TO SEE IF TWO 

SUCCESSIVE ITERATIONS ARE WITHIN LIMITS 

LHLD 

NEWT ♦ 2 


XCHG 



LHLD 

NEWT 


FSUB 

D»QLDT 




MOV 

A # D 


UP 

CQMPT 


FSUB 

D| 0D6DH 

I 8 * {OUDT «• NEWT) 

COMPTS FSUB 

D # 0D45H 

JTG, GET ABSOLUTE ERROR 

I 0,05 * DJFF 

JM 

NEXT 

I DIFF , GT , 0,05 

LHLD 

NEWR ♦ 2 


XCHG 

LHLD 

NEWR 


FSUB 

D.OLDR 


MOV 

A,D 


UP 

COMPR 


FSUB 

D» 0D6DH 

I FOR ABSOLUTE ERROR 

CQMPR? FSUB , 

D»0D71H 

UM 

NEXT 

;DIFF ',Gt', 1 FOOT 

jDIFFERENCE IS BETWEEN LIMITS 

AND HENCE ITERATIONS ARE COMPLETE 

ICQMPUTE 

DA TO DETERMINE 8F 


IDA * ARSINC <RP*COS<SP)*B#TF)/<RF»COS<SF>) ; 

LHLD 

SF + 2 


XCHG 

LHLD 

SF 

■■ . | 

1 

CCS 

fmul 

D, NEWR 


SHLD 

TEMP 


XCHG 



SHLD 

TEMP + 2 

iTEMP 5 NEWR * COS (SF ) 

LHLD 

0D07H 


XCHG 



LHLD 

0D05H 


COS 

FMUL 

D # ODOIH 


FMUL 

D.0D15H 


F MUL 

D# NEWT 


FDIV 

TEMP# D 


ARSIN 

FADD 

D.0D09H 

IBF * BP + DA 


5LIPIT 8F TO BE WITHIN 0 TO 360 DEGREES 



MOV 

A#D 



UK 

BFNEG 



FSUB 

QD4DH#D 

IBF - 6,28 


UP 

STORE 


b F M~ G J 

FADE 

0D4DH #D 


STORE I 

SHLD 

BF 



XCHG 




SHLD 

8F + 2 



j ##«***#***#*#***»#■##*##*##**##****#**#***#*•*##*#«#**#**#***** 



? INTRODUCING CORRECTIONS, 




{Confute tangent elevation Eg = g«tf«cos(Sf)/<2V) 
lhld sf + 2 

XCHG 

LHLD SF 

COS 

FMUL D| 0D55H 

FMUL DiNEWT 

FDIV 0D75H.D 

FDI V V# D 

SHLD EO 

XCHG 

SHLD EO * 2 

; »«»««•»«««»#«« » «#»«*»««« #«»«*«««»*«««««« «««««*«« ««««»«««*««««« 

{DRIFT CORRECTION 

{#««««««»#«««•»*««*««#«#«##**»*« «««»*«««« «*«««»«*»#««»*« MM** **** 

XCHG 

FMUL D,QD21H . • . , 

SHLD DRIFT {DRIFT = EO * DRIFT CONST, 

XCHG • 

SHLD DRIFT +2 


;# *««««« ««««## ##{*«»*»««»«« ««#«##««»««# ##««*« ##«#«««**« #mm**m*« * 

{MUZZLE VELOCITY AND BALLISTIC CORRECTIONS 

}##»*■&*■»*#»•»*■##'**•*■»*#**#«******»**■»* ##**«***««ihhhhh»**#mMM**M*< 


LHLD 

XCHG 

0D2BH 

IDELV 

LHLD 

0D29H 


FMUL 

0D31H t D 

{DELV * Cl 

SHLD 

XCHG 

TEMP 


SHLD 

TEMP + 2 


LHLD 

XCHG 

0D2FH 

{DELC 

lhld 

0D2DH 

,-DELC * C2 

FMUL 

0D35H # D 

FADD 

TiMP# P 


FaDD 

d.newt 


SHLD 

NEWT 


XCHG 


I TF s TF ♦ DELTM 

SHLD 

NEWT ♦ 2 

LHLD 

X CH 6 

0D2BH 

|DfeL V 



BNEG ? 
STOBJ 


LHLD 

0D29H 


FHUL 

0D39H # D 

5DELV * C3 

SHLD 

TEMP 


XCHG 

SHLD 

TEMP + 2 


LHLD 

GD2FH 


XCHG 


' 

LHLD 

0D2DH 


FHUL 

0D3DH, D 

|C4*BELC 

FADD 

D i TEMP 


fadd 

D|SF 


SHLD 

SF 


XCHG 

SHLD 

SF + 2 

1 SF ■ SF + DELVM 

! »»«« #»«####*■»•■&#***«###*»*■» 

l WIND 

CORRECTIONS, 



}#«*«»«««««**«*««#*»«»«#««#*«*#««#««««#«#«*«**««««*« **#mmm**m** 


LHLD 

SF + 2 

XCHG 


LHLD 

SF 

FADD 

D* EO 

COS 


FMUL 

D » WR 

FMUL 

DiNEWT 

FADD 

D,NEWR 

LHLD 

SF + 2 

XCHG 


LHLD 

SF 

COS 


SHLD 

COSSF 

XCHG 


SHLD 

COSSF + 2 

XCHG 


FMUL 

NEWR.D 

SHLD 

TEMP 

XCHG 


SHLD 

TEMP +2 

LHLD 

WL + 2 

XCHG 


LHLD 

WL 

FMUL 

DiNEWT 

fdiv 

TEMPiD 

FADD 

D,BF 

MOV 

A*D 

UM 

BNEG 

FSUB 

0D4DH # D 

UP 

STOB 

FADD 

QD4DH i D 

SHLD 

BF 

XCHG 


SHLD 

BF * 2 


iPHI = SF + EQ 

JWR *C0S(PHI)«TF 
?RF * RF + DELRW 


JRF * COS(SF) 


jMU # TF 

?WL«TF/<Rf*S0S(SF) 

JBF s BF + DELAW 

»BF TO pE 0. TO 360 DEGREES 
J BF - 6,28 



COSSF + 2 


LHLD 
XCHG 
LHLD COSSF 

FMUL D,N£WT 

fmul D,WR 

FD I V V • D 

FADD d,newt 

SHLD NEWT 

XCHG 

SHLD NEWT + 2 

LHLD SF + 2 

XCHG 

LHLD SF 

FADD D » E 0 

SHLD TEMP 

XCHG 

SHLD TEMP + 2 

XCHG 


Cos 

SHLD 

COSP 

XCHG 

SHLD 

COSP + 2 

LHLD 

TEMP + 2 

XCHG 

LHLD 

TEMP 

SIN 

FMUL 

D # COSP 

FMUL 

D# V 

SHLD 

TEMP 

XCHG 

SHLD 

TEMP + 2 

LHLD 

COSSF + 2 

XCHG 

LHLD 

COSSF 

FMUL 

WR.D 

FDIV 

temp,d 

FSUB 

D , SF 

SHLD 

SF 

XCHG 

SHLD 

SF + 2 

XCHG 

FADD 

EQ# D 

SHLD . 

GUNEL 

XCHG 

SHLD 

GUNEL + 2 

JMP 

TIMER 


iDELTW=WR*TF*CQS<SFj/V 
JTF s IF + DELTW 

? COS (PHI ) 

# S I M < PH | ? 

J COS (PHI )*SIN(PHI ) 

■ V 

;SF * SF • DELVW 

JRETURNS TO MAIN ROUTINE 


|WHEN ITERATIONS TO CONTINUE TO GET SUCCESSIVE VALUES OF RANGE 
JAND TIME OF FLIGHT WITHIN TOLERANCES 


NEXT I f BINS NpNR* OLDR 
DCR A 

UM LOST 

-•TRANS NEWT.OLDT 
TRANS NEWRiQLDR 


IOLDR s NEWR 

J DOES NOT CONVERGE IN *5 
I ITERATIONS 

?oldt s newt 

IOLDR = NEWR 


QD03M 


LHLD 
XCHG 
LHLD ODOlH 


^UL D.ODliH 

FMUL D#NEWT 

niv NEWR* D 

A RS I N 

SHLD DV 

XCHG 

SHLD DV ♦ 2 

XCHG 

FADD D , 0D05H 

SHLD SF 

XCHG 

SHLD SF + 2 

RIMIT SF TO BE WITHIN 10 TO 

XCHG 

MOV A, D 

JM SFNEG 

FSUB D,0D5lH 

CM SFH I 

FSUB D.0D51H 

SHLD SF 

XCHG 

SHLD SF + 2 . • 

?SF ADJUSTED BETWEEN 10 TO 6 

LHLD NEWR + 2 

XCHG 

LHLD NEWR 

FSUB D * QD8DH 

JM LOOP 

LHLD SF + 2 

XCHG 

LHLD SF 

FSUB D * 9D59H 

JP SFTEN 

FADD D » QD69H 

JM SFETY 

JMP LOOP 

SFTENS TRANS 0D59H,SF 
JMP LOOP 

SFETY 5 TRANS 0D5DH,SF 

JMP LOOP 

SFnEg: FaDD d,dv 

, JMP LOOP* 

SFH I S TRANS QD61W,SF 

JMP LOOP " ' 

* WHEN ROOT RETURNED IS ZERO 

L00P2I LHLD SF * 2 
XCHG 


JRP * S 


» SF * SP + DV 
SO DEGREES 

ilF SF # LT', C 
I 90 - SF 
JSF * G T , 90 
{RESTORE SF 

deg, when Range below 3000 0 feet 

i 30000 * RF 

t 10 - SF 

i 70 + 10 » SF = 80 f SF 

! SF s 10 

J SF = 80 

; SF s SF + DV ^ DV 

i SF a 89 

TRY AGAIN BY ALTERING SF 



lhld 

FftDC 

F3UB 

JK 

JKP 

L i" ■ S T I X R A 
NOV 

kcv 

SHLD 

XChG 

SHLD 

JKP 


SF 

D# 0D6SH 
D * 0D51H 
LOST 
LOOP 
A 

H# A 
L » A 
NEWR 

NEWR + 2 
TIMER 


I SF .* SF + 5 
190 - SF 


I ROUT I NE TO DETERMINE SQUARE ROOT BY NEWTON * RAPHSON METHOD 

i********************************************************.***, 

inr'i^r^^T’ DEPENDS On WHETHER RP/V IS Ge'. V/(G*SIN£SF) ) 
?OR IT IS » LT , V/ ( G#S IN (SF ) ), ’ * ’ V/ ' b blN(SF)) 

,RP/V ,GE * V/(G#SIN(SF ) ) IMPLIES ROOT IS<-B/(2*A> ) + SQROOT 
IRP/V .IT*. V/(G*SIN(8F)) IMPLIES ROOT IS <-B/(2#A>> - SQROOT 


ROOTS? XRA A 


STA 

FLAG 

LHLD 

BB + 2 

XCHG 


LHLD 

BB 

FMUL 

D#B8 

SHLD 

TEMP 

XCHG 


SHLD 

TEMP ♦ 2 

LHLD 

AA + 2 

XCHG 


LHLD 

AA 

FMUL 

D,0DB1H 

FMUL 

D,CC 

FSUB 

D, TEMP 

SHLD 

PP 

XCHG 


SHLD 

PP + 2 

LHLD 

QDQ3H 

XCHG 


Lhld 

0D01H 

FDIV 

V»D 

SHLD 

TEMP 

XCHG 


SHLD 

TEMP | 2 

lhld 

A A ♦ 2 

XCHG 



jdecides add or sub sgroot 


*4 * AA 
JJ « AA « CO 
?BB«BB - 4#AA#CC 


$ CHECK RP/V DIFF V/G*SIN(SF 



Fd X v 
Fsue 
jp 

fslb 


D * V 

D'TEpp 

X1R 

B » 0D6DH 


JHP 

XiR: HVI 

IB CPs SHLD 
XCHG 
SHLD 
XCHG 

fmul 

FSUB 

FSUB 

SHLD 

XCHG 

SHLD 

LHLD 

XCHG 

LHLD 

FMUL 

FD IV 

SHLD 

XCHG 

SHLD 

FSUB 

HO V 

UP 

FSUB 
A USD! FSUB 
JP 

Trans 

JMP 

QUTLi LDA 
JNZ 
LHLD 
XCHG 
LHLD 
FADD 
FDiV 
FD I V 
FSUB 
SHLD 
XCHG 
SHLD 
, JMP 

X1VALI LHLD 
XCHG 
LHLD 
FSUB 
FADD 
F0IV 
FDIV 


ILOP 

A*1 

X 

X + 2 


D,X 
PPiD 
D * X 
TEMP 

TEMP +2 
X + 2 

X 

D» 0D75H 
D i TEMP 
XXI 

XXI + 2 
X i D 
A* D 
ABSD 
D 1 OD 6 DH 
D, OD9DH 
OUTL 
XX1,X 
ILOP 

flag 

xival 

B B + 2 
BB 

XXI, D 
A,D 

0D75W t D 
D.0D6DH 
NEWT 

NEWT 4 2 
ROQTR 
88'* 2 

BB 

D # O 06 DH 

D.XXl 

A,D 

nD^Bu *n 


JRP/V * V/G#SIN<SF) 


iRP/V ,LT , V/(G#SIN(SF ) ) 
WHENCE INITIAL GUESS IS 
JV/( G#SJN(SF) ) - Rp/V 


'FLAG SET’, ADD SQRQOT 
iX = INITIAL GUESS 


JX*X - PP = F ( X ) 

;x * F(X) 


} 2 * X a F» (X) 
*X*F(X)/F»(X) 


» (XXI - X) a diff 

'FOR A BS , V AL UE 

ITERATIONS COMPLETE 
s XXI 

fADD SQROOT 

J SUBTRACT SQROOT 
18 + SQRQOT 

» (B+SqROOT )/( 2 *A ) 

,*0 (B+SQROOT > / ( 2 *A ) ) 


id ■ 8 s 


ft J ... W M M M #>1 a W 1 A * mi . A . 



f% V ! "1 W 

ariLL 

NEWT + 2 


JMP 

RQCTR 


£ i\ D 01 

PROGRAM 


OS 

1 

«. T a : 

fT C 
fc»i *><* 

i 

r J 

cs 

1 

K j 

I" 8 

1 

l t 

CS . 

1 


ns 

1 

' J. 5 

D3 

1 

2 t 

DS 

1 

h i\ : 

CS 

1 

r L fi G s 

E S 

1 

hr 1 ! 

DS 

• 2 

S P c 5 

DS 

2 

S r 3 ! 

DS 

2 

S F 4 ! 

DS 

2 

IPT1: 

DS 

2 

IPT2! 

DS 

2 

IPT3! 

DS 

2 

1PT4: 

DS 

2 

STACl ! 

DS 

2 

i£ T A C 2 S 

DS 

2 

ST ACS J 

DS 

2 

STAC A 5 

DS 

2 

CLOCK* 

DS 

4 

WR j 

DS 

4 

bit 

DS 

4 

V! 

DS 

4 

CLCRl 

DS 

4 

SF S 

DS 

4 

CLDTS 

DS 

4 

AA j 

DS 

4 

BB S 

DS 

4 

CCI 

DS 

4 

T BhP s 

DS 

4 

F P : 

DS 

4 

yt 

DS 

4 

xx a * 

DS 

4 

N E K T 1 

DS 

4 

KEKRI 

DS 

4 

BF 1 

DS 

4 

£0 * 

DS 

4 

DRIFT* 

DS 

4 

CCSSF S 

DS 

4 

CCSP J 

DS 

4 

CUR EL! 

DS 

4 

DV * 

DS 

4 

LR*S 

DS 

4 

L;S X ! 

DS 

4 

DBX* 

DS 

4 

CL CK 1 1 

DS 

4 

ACPI! * 

DS 

4 



M »„i h ,X * 

0 

r j « 

D S 

a ■ l. ;t s 

OS 

» 1 1-.' ,i i 

cs 

' r « 

l s 

A *' i * V # 

cs 

A l - C, 4 « 

cs 

^ J f * A’ # 

L S 

z 1 'i £ 2 

DS 

/ ;i £ < 

c s 

. h 3 * 

cs 

:!i 3 ' * 3 ; 

r c 

*W J5 

A w '2 • 

> ! ^ < » V f 

L S 

a i B j : 

DS 

A , t c : 

c s 

a 1 u o ; 

DS 

A A ft { 

CS 

>. A £ 1 

DS 

a A & ; 

CS 

E B ft s 

CS 

8 B E ! 

DS 

EBB: 

DS 


*4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 
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This book is to be returned on the 
date fast stamped. 
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